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THE PUPIL DETECTION METHOD AND SHAPE DESCRIPTOR EXTRACTION 
METHOD FOR A IRIS RECOGNITION, IRIS FEATURE EXTRACTION 
APPARATUS AND METHOD, AND IRIS RECOGNITION SYSTEM AND 

METHOD USING ITS 

5 

Technical Field 

The present invention relates to a biometric 
technology based on a pattern* recognition and a image 

10 processing; and, more particularly, to a pupil detection 
method and a shape descriptor extraction method for an iris 
recognition that can provide a personal identification 
based on ari. iris of . an eye, .an . iris .feature extraction 
apparatus, and method , and\. iris recognition system and method 

15 using the same, and a computer-readable recording • medium 
that records programs impliementing the methods. 

Packqround hx\. 

20 Conventional methods for identifying a person, e.g., a 

password and a personal identification number, cannot 
provide accurate and reliable personal identification in an 
information society that is getting highly developed, due 
to stealth or lost of the password and the identification 

25 number, and cause side effects according to a reverse 
function. 

Particularly, it is predictable that rapid development 
of an internet environment and increase of tthe electronic 
commercial cause enormous mental blow and material damage 
30 to a personal or an organization using only those 
conventional identification method. 

Since among various biometric methods, the iris is 
broadly known as most effective in a view of identity, 
invariance and stability, and a failure rate of the 
35 recognition is very low, the iris is applied to a field 
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that requires hi^h security. , 

Generally^ in method for identifying a person using 
the iris, it is indispensable to detect speedily a- pupil 
and the iris for real-time iris recognition from an image 
5 signal of an eye of the person. 

Hereinafter, features of the* iris and a conventional 
method for the iris recognition will be described. 

In a process for precisely dividing the pupil from the 
iris by detecting a pupil boundary, it is very important to 
10 achieve a feature point and a normalized feature quantity 
regardless of a pupillary dilation without allocating the 
same part of the iris analysis region to the same 
coordinates when the image is analyzed. 

Also, the feature point of the iris analysis region 
15 reflects an iris fiber, a structure of layers and a 
defection of a connection state. Because the structure 
. affects to a function and reflects intei^rity, the structure 
indicates a resistance of an organic and a genetic factor. 
As related signs, there are lacunae^ crypts, defect signs 
20 and rarifition and so on. 

The pupil is located in the middle of the iris and 
iris collarette that is an iris frill having a sawtooth 
shape, i.e., autonomic nerve wreath in the iridology, is 
located at l'-2 mm distance from a pulillary. Inside of the 
25 collarette is an annuls iridis minor and outside of the 
collarette is an annuls iridis major. The annulus iridis 
major includes iris furrows that are a ring-shape 
prominence concentric to the pulillary. The iris furrows 
are referred to as a nerve ring in the iridology. 
30 In order to. use an iris pattern based on a clinical- 

experience of the iridology as the feature point, the iris 
analysis region is divided into 13 sectors and each sector 
is subdivided into 4 circular regions based on a center of 
, the pupil. 

35 The iris recognition system extracts an image signal 
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from the iris, transforms the image signal to specialized 
iris data, searches identical data to the specialized iris 
-data in a database and compares the searched data to the 
^ specialized iris data, to thereby identify the person for 
5 acceptance or refusal. 

It is important to search a statistical texture, i.e.,' 
an iris, shape, in the iris recognition system. Features 
that a person recognizes the texture are periodicity, 
directionality and randomness in a cognitive science. 
10 Statistical feature of the iris includes the degree of 
freedom and sufficient identity to identify a person. An 
individual can be identified based on the statistical 
feature. 

Generally, in the conventional pupil extraction method 

15 of the conventional iris recognition system proposed by 
Daugman, a circular projection is obtained at every 
location of the image and a differential value of the 
circular projection is calculated, and then the largest 
value obtained by calculating the differential value based 

20 on Gaussian convolution is estimated as the boundary. Then, 
a location that the circular boundary component is the 
strongest is obtained based on the estimated boundary, to 
thereby extract the pupil from the iris image. 

However, it takes long time to extract the pupil 

25 because the projection for the whole image and the 
differential calculation increase operation numbers. 
Because it assumed that there is a circular component, it 
cannot be detected that there is no circular component in 
the conventional method. 

30 Also, the pupil detection must be processed before the 

iris recognition, and fast pupil extraction is required for 
real-time iris recognition. . However, if a light source 
exists in the pupil, ^n inaccurate pupil boundary is 
detected due to infrared rays. Because the above problem, 

35 the iris analysis region must be whole image except light 

3 
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origin Region. Therefore, the accuracy of the analysis is 
decreased. 

In particular, a method for dividing a frequency 
region based on a filter bank and extracting the 
5 statistical feature is generally used in the iris feature 
extraction. Gabor filter or Wavelet, filter is used. The 
Gabor filter can divide the frequency region effectively, 
arid the Wavelet filter can divide the frequency region on 
consideration of a human eyesight feature. However, 

10 because the above methods require many operations, i.e., it. 
needs much time, the above methods are not appropriate for 
the iris recognition system. In detjail, because much time 
• and cost are needed to develop the iris recognition system 
and the recognition operation cannot be operated rapidly, 

15 the method for extracting the statistical feature is not 
effective. Also, because the feature value is .not 
rotation-invariant or scale-invariant, there is a 
limitation that the feature value is rotated and compared 
in order to search the converted texture. 

20 ' However, in the case of the shape, it is possible to 

search the boundary by expressing in direction, and to 
express and search the shape of the image regardless of 
change, motion, rotation and scale of the shape by using 
various transformations. Therefore, it is desirable to 

25 preserve an iris boundary shape or an efficient feature of 
a part of the iris. 

A shape descriptor is based on a lower abstractic^ 
level description that can be automatically extracted, and 
is a basic descriptor that human can indicate from the 

30 image. There are two well-known shape descriptors adopted 
by experiment model (XM) that is a standard of Motion 
Picture' Expert Group-7 (MPEG-7) . The first shape 

descriptor is Zernike moment shape descriptor. A Zernike 
basis function is prepared in order to get distribution of 

35 various shapes in the image and the image having a 
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.predetermined sized is ^projected to the basis function, and 
the projected value is used as the Zernike moment shape 
descriptor. The second shape descriptor is Curvature scale 
§pace descriptor. A low frequency-pass filtering of the 
5 contour extracted from the image is' performed, a change of 
inflection point existing on the contour is expressed in a 
scale space, and a peak value and the location of the 
inflection point are expressed as a two-dimensional vector. 
The two-dimensional vector is used as a Curvature scale 

10 space descriptor. 

Also, according to an image matching method using the 
conventional shape descriptor, a precise object must be 
extracted from the image in order to search a model image 
having the shape descriptor similar with the shape 

15 descriptor, of a query image. Therefore, it is a drawback 
that the model image cannot be searched if the object is 
not extracted precisely. 

.Therefore, it is required that a method for developing 
similar group database indexed based on a similarity shape 

20 descriptor, e.g., the Zernike moment shape descriptor or 
the Curvature scale space shape descriptor, and searching 
an indexed iris group having similar shape, descriptor with 
the query image from the database. In particular, the 
above method is very effective to 1:N identification (N is 

25 a natural number) . 

Disclosure ' . 

Technical problem .of the Invention 

30 It is, therefore, an object of the present invention 

to provide a method for extracting a pupil in real time and 
an iris .feature extraction" apparatus using the same for the 
iris recognition that, is not sensitive to illumination 
lighting to an eye and has high accuracy, and a computer- 

35 readable recording medium recording a program that 



wo 2005/024708 



PCT/KR2004/002285 



implements the methods. 

Also^ it is another object of the present invention to 
provide a method for extracting a shape descriptor which is 
invariant to motion, scale, illumination and rotation, a 
5 method for developing a similar groujJ database indexed by 
using the shape descriptor . and searching the index iris* 
group having a similar shape descriptor with the query 
image from the database, and an iris feature extracting 
apparatus using the same, an iris recognition system and a 

10 method thereof, and a computer-readable recording medium 
recording a program that implements the methods . 

Also, it is still another object of the present 
invention to provide a method for developing an iris shape 
database according to a ..dissimilar shape descriptor by 

15* measuring dissimilarity of a similar iris shape group 
indexed by the shape descriptor extracted by a linear shape 
descriptor extraction method and searching the index iris 
group having the shape descriptor matched to the query 
image from the database, and an iris feature extracting 

20 apparatus using the same, an iris recognition system and a 
method thereof, and a computer-readable recording medium 
recording a program that implements the methods.* 

Other objects and benefits of the present invention 
will be described hereinafter, and ^will be recognized 

25 according to an embodiment of the present invention. Also, 
the objects and the benefits of the present invention can 
be implemented in accordance with means and ^combinations 
shown in claims of the present invention. 

30 Technical solution of the Invention 

In accordance with an aspect of the present invention, 
there is provided a method for detecting a pupil for iris 
recognition, including the steps of: a) detecting light 
sources in the pupil from an eye image as two reference 

35 points; b) determining first boundary candidate points 



wo 2005/024708 PCT/KR2004/00228S 

located between the iris and the pupil of the eye image^ 
which cross over a straight line between the two reference 
points; c) determining secofid boundary candidate points 
located between the iris and the pupil of the eye image, 
5 which cross over a perpendicular bisector of a • straight 
line . between the first boundary candidate points; and" d) 
determining a location and a size of the pupil by obtaining 
a radius of a circle and coordinates of a center of the 
circle based on a center candidate point, wherein' the 

10 center candidate point is a center point of perpendicular 
bisectors of straight line between the neighbor boundary 
candidate points, to thereby detect the pupil. 

In accordance with another aspect of the present 
invention, there is provided a method for extracting a 

15 shape descriptor for iris recognition, the method including 
the steps of: a) extracting features of ah iris under a 
scale-space and/or a scale illumination; b) normalizing a 
low-iDrder moment with a mean size and/or a mean 
illumination, to thereby generate a Zernike moment which is 

20 size-invariant and/or illumination-invariant, based on the 
low-order moment; and c) extracting a shape descriptor 
which is rotation-invariant, size-invariant and/or 
illumination-invariant, based on the Zernike moment. 

The above method further includes the steps of: 

25 establishing an indexed iris shape grouping database based 
on the shape descriptor; and retrieving an indexed iris 
shape group based on an iris shape descriptor similar to 
■ that of a query image from the indexied iris shape grouping 
database. 

30 In accordance with another aspect of the present 

invention, there is provided a method for extracting a 
shape descriptor for iris recognition, the method including 
the steps of: a) extracting a skeleton from the iris; b) 
thinning the skeleton, extracting straight lines by 

35 connecting pixels in the skeleton, obtaining a line list; 

7 
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and c) normalizing the line list and setting the normalized , 

^ , , . 

line, list as a shape descriptor. 

The above method furt;her includes the steps of: 

establishing an iris, shape database of dissimilar shape 

5 descriptor by measuring dissimilarity of the images in an 

indexed similar iris shape group based on the shape 

descriptor; and retrieving an iris shape matched to a query 

.image from the iris shape database. 

In accordance with another aspect of the present 

10 invention,/ there is provided an apparatus for extracting a 
feature of iris, including: image capturing unit ■ for 
digitalizing and quantizing an image and obtaining an 
appropriate image for iris recognition; reference point 
detecting unit for detecting reference points in a pupil 

15 from the image, and detecting an actual center point of the 
pupil; boundary detecting unit for detecting an inner 
boundary between the pupil and the iris and an outer 
boundary between the iris and a sclera, to thereby extract 
an iris image from the image; image coordinates converting 

20* unit for converting a coordinates of the iris image from a 
Cartesian coordinates system to a polar coordinates system, 
and defining the center point of the pupil as an origin 
point of the polar coordinates system; image analysis 
region defining unit for classifying analysis regions of 

25 the iris image in order to use an iris pattern as a feature 
point based on clinical experiences of the iridology; image 
smoothing unit* for smoothing the image by performing a 
scale space filtering of the analysis region of the iris 
image in order to clearly distinguish a brightness 

30 distribution difference • between neighboring pixels of the ' 
image; image normalizing unit for normalizing a low-order 
moment used for the smoothen image as a mean size; and 
shape descriptor extracting unit for generating a Zernike 
moment based on - the feature point extracted in a scale 

35 space and a scale illumination, and extracting a shape 
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descriptor which is. rotation-invariant and noise-resistant 
by using Zernike moment. 

.The above apparatus further includes reference value 
storing unit for storing a reference value as a template by 
5 comparing a stability of the Zernike moment and a 
similarity of Euclid distance. 

In . accordance with another . aspect of the present 
invention, there is provided a system for recognizing an 
iris, including: image capturing unit for digitalizing and 

10 quantizing an image and obtaining an appropriate image for 
iris recognition; reference point • detecting unit for 
. detecting reference points in a pupil from the image, and 
detecting an actual center point of the pupil; boundary 
detecting unit for detecting an inner boundary between the 

15 pupil and the iris and an outer boundary between the iris 
and a sclera, to thereby extract an iris image from the 
image; image coordinates converting unit for converting a 
coordinates of the iris image from a Cartesian coordinates 
system to a polar coordinates system, and defining the 

20 center point of the pupil as an origin point of the polar 
coordinates system; image analysis region defining unit for 
classifying analysis regions of the iris image in order to 
use an iris pattern as a feature point based on clinical 
experiences of the iridology; image smoothing unit for 

25 smoothing the image by performing a scale space filtering 
of the analysis region of the iris image . in order to 
clearly distinguish a brightness distribution difference 
between neighboring pixels of. the image; image ' normalizing 
unit for normalizing a low-order moment used for the 

30 smoothen image as a mean size; shape descriptor extracting 
unit for generating a Zernike moment based on the feature 
point extracted in a. scale ^ space and a scale illumination; 
and extracting a* shape descriptor which is* rotation- 
invariant and noise-resistant by using Zernike moment; 

35 reference value storing unit for storing a reference value 
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as a template hy comparing a stability of the ^Zernike 

moment and a similarity . of Euclid distance; and 

verifying/authenticating unit for verifying/authenticating 

the iris by matching the feature quantities between models 

5 each of which represent the stability and the similarity of 

the Zernike moment of the query iris image in statistical. 

In accordance with another aspect of the present 

invention/ there is provided a method for extracting a 
■ .-'IP 
feature of an iris, including the steps pf : a) digitalizing 

10 and quantizing an image and obtaining an appropriate image 
for iris recognition; b) detecting reference points in a 
pupil from" the image, and detecting an actual center point 
of the pupil; c) detecting an inner boundary between the 
pupil and the iris and an outer boundary between the iris 

15 .and a sclera, to thereby extract an iris image from the 
image; d) converting a coordinates of the iris image from a 
Cartesian coordinates system to a polar coordinates ' system, 
and defining the center point of the pupil as an origin 
point of the polar coordinates system; e) classifying 

20 analysis regions of the iris image in order to use an iris 
pattern as a feature point based on clinical experiences of 
the iridology; f) smoothing the image by performing a scale 
space filtering of the analysis region of the iris image in 
order to clearly distinguish a brightness distribution 

25 difference between, neighboring pixels of the image; g) 
normalizing the image by normalizing a low-order moment 
with a mean size, wherein the low-order moment is used for 
the smoothen image; and h) generating a Zernike moment 
based on the feature point extracted in a scale space and a 

30 scale illiamination, and extracting a shape descriptor .which 
is rotation-invariant and noise-resistant by using Zernike 
moment. 

The above method further includes the step of i) 
storing a reference value as a template by comparing a 
35 stability of the Zernike moment and a similarity of Euclid 

10 
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. distance. . : 

•In . accordance /with another aspect of the' present, 
invention, there is provided a method for recognizing an 
iris, * including the steps of:' a) digitalizing and 

5 quantizing an image and obtaining an .appropriate image for 
iris recognition; b) detecting reference points in a pupil 
from the image, and detecting an actual center point of the 
. pupil; c) detecting an inner boundary between .the pupil and 
the iris and an outer boundary between the iris and a 

10 • sclera, to thereby extract an iris image from the image; d) 
converting a coordinates of the iris image from a Cartesian 
coordinates system to a polar coordinates system, and 
defining the center point of the pupil as an origin point 
of the polar coordinates system; e) classifying analysis 

15 regions of the iris image in order to use an iris pattern 
as a feature point based on clinical experiences of the 
iridology; f) smoothing the image by performing a scale 
space filtering of the analysis region of the iris image in 
order to clearly distinguish a brightness distribution 

20 difference between neighboring pixels of the image; g) 
normalizing the image by normalizing a low-order moment 
with a mean size, wherein the low-order moment is used for 
the smoothen image; h) generating a Zernike moment based on 
the feature point extracted in a scale space and a scale 

25 illumination, and extracting a shape descriptor which is 
rotation-invariant and noise-resistant by using Zernike 
moment; i) storing a reference value as a template by 
; comparing . a stability of the Zernike moment and a 
similarity of Euclid distance; and j) 

30 verifying/authenticating the iris by matching the feature 
quantities between models each of which represent the 
stability and the similarity of the Zernike moment of the 
query iris image in statistical. 

In accordance with another aspect of the present 

35 invention, there is provided a computer readable recording 
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medium storing program for executing a method for detecting 
a pupil for \ iris recognition, the method including the 
steps of: a), detecting iight sources in the pupil from an 
eye image as two reference points; b) determining first 
5 boundary candidate points located between the iris and the 
pupil of the eye image, which . cross over a straight line 
between the two reference points; c) determining second 
boundary candidate points located between the iris and the 
•pupil, of the eye image, which cross over a perpendicular 

10 bisector of a straight line between the first boundary 
candidate points; and d) determining a location and a size 
of the pupil by . obtaining a radius of a circle and 
coordinates of a center of the circle based on a center 
candidate point, wherein the center candidate point is a 

15 center point of perpendicular bisectors of straight line 
between the neighbor boundary candidate . points, to thereby 
detect the pupil. 

In accordance with another aspect of the present 
invention, there is provided a computer readable recording 

20 medium storing program for executing a method for 
extracting a shape descriptor for iris recognition, the 
method including the steps of: a) extracting a feature of 
iris under a scale-space and/or a scale illumination; b) 
normalizing a . low-order moment with a mean size and/or a 

25 mean illumination, to thereby generate a Zernike moment 
which is size-invariant and/or illumination-invariant, 
based on the low-order moment; and c) extracting a shape 
descriptor which is rotation-invariant, size-invariant 
and/or illumination-invariant, based on the Zernike moment. 

30 The above computer readable recording medium further . 

includes the steps of: establishing an indexed iris shape 
grouping database based on the shape descriptor; and 
retrieving an indexed iris shape group based on an iris 
shape descriptor similar to that of a query image from the 

35 indexed iris shape grouping database. 

12 
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In accordance with another aspect *o^ th.e. pres^ent 
invention, there is provided a computer readable recording 
m^ium storing program for executing a method for 
extracting a shape descriptor for iris recognition, the 

5 method including the steps of: a) extracting a skeleton 
from the iris; b) thinning the skeleton, extracting 
straight lines by connecting pixels in the skeleton, 
obtaining a line list; and c) normalizing the line list and 
setting the normalized line list as a shape descriptor, 

10 The above computer readable recording medium further 

includes the steps of: establishing ah iris shape database 
of dissimilar shape descriptor by measuring dissimilarity 
of the images in an indexed . similar' iris, shape group based 
on the shape descriptor; and retrieving an iris shape 

15 matched to a query image from the iris shape database. 

In accordance with another aspect of the present 
inviention, there is. provided a computer readable recording 
medium - storing program for executing a method for 
extracting a feature of iris, the method including the 

20 steps of:, a) digitalizing and quantizing an image and 
obtaining an appropriate image for iris recognition; b) 
detecting reference points in a pupil from the image, and 
' detecting an actual center point of the pupil; c) detecting 
an inner boundary between the pupil and the iris and an 

25 outer boundary between the iris and a sclera, to thereby 
extract an iris image from the image; d) converting a 
coordinates of the iris image from a Cartesian coordinates 
system to' a polar coordinates system, and defining the. 
center point of the pupil as an origin point of the polar 

30 coordinates system; e) classifying analysis regions of the 
iris image in order to use an iris pattern as a feature 
* point based on clinical experiences of the iridology; f) 
smoothing the image by performing a scale space filtering 
of the analysis region of the iris image in order to 

35 clearly distinguish a brightness distribution difference 
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between neighboring pixels of the image; g) normalizing the 
image by normalizing a ' low-order moment with^a mean size, 
whjerein the low-order moment is used for the smoothen 
image; and h) generating a Zernike moment based on the 
5 feature point extracted in a scale space and a scale 
illumination, and extracting a shape descriptor which is 
rotation-invariant and noise-resistant by using Zernike 
moment . 

The above computer readable recording medixam further 

10 includes the step of: i) storing a reference value as a 
template by comparing a stability of the Zernike moment and 
a similarity of Euclid distance. 

In accordance with another aspect of the . present 
invention, there is provided a computer readable recording 

15 medium which is recorded program for executing a method for 
recognizing an iris, the method including the steps of: a) 
digitalizing and quantizing an image and obtaining an 
appropriate image for iris recognition; b) detecting 
reference points in a pupil from the image, and detecting 

20 an actual center point of the pupil; c) detecting an inner 
boundary between the pupil and the iris and an outer 
boundary between the iris and a sclera, to thereby extract 
an iris image from the image; d) converting a coordinates 
of the iris image from a Cartesian coordinates system to a 

25 polar coordinates system, and defining the center point of 
the pupil as an origin point of the polar coordinates 
system; e) classifying analysis regions of the iris image 
in order to use an iris pattern as a feature point based on 
clinical experiences of the iridology; f) smoothing the 

30 image by performing a scale space filtering of the analysis 
region of the iris image in order to clearly distinguish a 
brightness distribution difference between neighboring 
pixels of the image; g) normalizing the image by 
normalizing a low-order moment with a mean size, wherein 

•35 the low-order moment is used for the smoothen image; h) 
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. generating "a^ Zernike moment based on the feature point 
extracted in a scale space and a scale illumination, and 
extracting a shape descriptor which is rotation-invariant 
and noise-resistant by using Zernike moment; i) storing a 

5 • reference value as a template by comparing a stability of 
the Zernike moment and a. similarity of Euclid distance; and 
j) verifying/authenticating the iris by matching the* 
feature quantities between models each of which represent 
tihe stability and the similarity of the Zernike moment of 

10 the query iris image in statistical. 

The present invention provides an identification 
system which identifies a person or discriminate the person 
from others based on the iris of an eye quickly and 
precisely. The identification system acquires an iris 

15- pattern image for iris recognition, detects an iris and a 
pupil quickly for real-time iris recognition, extracts the 
unique features of the iris pattern by solving the problems 
of a non-contact iris recognition method, i.e., variation 
in the image size, tilting and moving, and utilizes the 

20 Zernike moment having the visional recognition ability of a 
human being, regardless of motion, scale, illumination and 
rotation. 

For the identification system, the present invention 
acquires an image appropriate for the iris recognition by 

25 computing the brightness of an eyelid area and the pupil 
location based on the iris pattern image, performs 
diffusion filtiering in order to remove noise in the edge 
area of an iris pattern image obtained by carrying out 
Gaussian blurring, and detects the pupil in real-time more 

30 quickly by using a repeated threshold value changing method. 
Since pupils have a different curvature, their radiuses are 
obtained by using a Magnified Greatest Coefficient method. 
Also, the central coordinates of a pupil is obtained by 
using a bisection method and then the distance from the 

35 center of the pupil to the radius of the pupil is obtained 
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in the counter clock-wise. Subsequently^ the precise 
boundary is. detected by taking the x-axis as a rotational 
angle and the y-axis* as a distance from the center to the 
radius of the pupil and expressing the result in a graph. 

5 Also, the iris features are" extracted through a scale- 

space filtering. Then, the Zernike moment having an 
invariant feature is generated' by using a low-order moment 
and the low-order moment is normalized with a me^n size in * 
order "to obtain features that are not chang.ed by the size, . 

10 illumination and rotation. The Zernike moment is stored as 
a reference value. The identification system 

recognizes/identifies an object in the input image through 
a feature quantity matching between models reflecting the 
similarity of the reference value, the stability of the 

15 Zernike moment of the input image, and the feature quantity 
in probabilities. Herein, the identification system can 
identify an iris of a living person quickly and clearly by 
combining the Least Square (LS) and Least Media of Square 
(Lmed) algorithms. 

20 To be more specific, the present invention directly 

acquires a digitalized eye image by using a digital camera 
instead of using a general video camera for identification, 
selects an eye image appropriate for recognition, detects a 
reference point within a pupil, defining a boundary between 

25 the iris and the pupil of the eye, and then defining 
another circular boundary between the iris and a sclera by 
using an arc that does not necessarily form a concentric 
circle with the pupil boundary. In other words, the 
identification system directly acquires a digitalized eye 

30 image by using a digital camera instead of a general video 
camera for identification, selects an eye image appropriate 
for recognition, detects a reference point within . the pupil, 
detecting a pupil boundary between the iris and the pupil 
of the eye, detecting a pupil region by acquiring the 

35 center coordinates and the radius of the circle and 
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' determining the location and size* of the pupil/ detects 
an external area between the iris region and the sclera 
region' by using an^ arc ^that does not ne^cessarily form a 
concentric circle With the pupil boundary, 

5 A pplar coordinate system is established and the 

center of the circular pupil boundary of the iris pattern 
image is put in the origin of the polar coordinate system. 
Then, an annular analysis region is defined within the iris. 
The analysis region* appropriate for recognition does not 

10 include pre-selected parts , e.g., the eyelid, the eyelashes 
or a part can be blocked off by mirror reflection from 
illumination. The iris pattern image in the analysis 
region is transformed into a polar coordinate system and 
goes through 1-order scale-space filtering that provides 

15 the same pattern regardless of the size of the iris pattern 
image by using a Gaussian kernel with respect to a one- 
dimensional iris pattern image of the same radiuses around 
the pupil. Then, an edge, which is a zero-crossing point, 
is obtained,, and the iris features are extracted in two- 

20 dimensional by accumulating the edge by using an overlapped 
convolution window. This way, the size of data can be 
reduced during the generation of an iris code. Also, the 
extracted iris features can make a size-invariant Zernike 
moment which is rotation-invariant but sensitive to size 

25 and illumination as normalizing the moment into a mean size 
by using the low-order moment in order to obtain a feature 
quantity. If a change in a local illumination is modeled 
into a scale illumination change and the moment is 
normalized into a mean brightness, an illumination- 

30 invariant Zernike moment can be generated. . A Zernike 
moment is generated based on the feature point extracted 
from the scale space and scale illumination and stored as a 
reference value. At a recognition part, an object in the 
iris image is identified by matching the feature quantity 

35 between models reflecting the reference value, stability of 
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the Zernike moment * and similarity between feature 
quantities in probability. Wherein, the iris recognition 
4.S Verified by combining the LS and 'the Lmeds methods.. * • ^ 
In accordance with the present invention, the feature 

5 quantity that is invariant to a local illumination change 
is generated by changing a local Zernike moment based on 
biological facts that a person focuses at the main feature 
point when a person recognize the object. Therefore, an 
image of the eye must be' acquired as a digital form 

10 appropriate to analyze. Then, an iris region of the image 
is defined and separated. The defined region of the iris 
image is analyzed, and to thereby generate the iris feature. 
A moment based on the feature generated for a specific iris 
is generated and stored as a reference value- In order to 

15 obtain an outlier, the moment of the input image is 
filtered using the similarity and the stability used for 
probability object recognition and thign 'is matched to the 
stored reference moment. The outlier allows the system to 
confirm or disconfirm the identification of the person and 

20 evaluate confirm level of the decision. Also, a 

recognition rate can be obtained by discriminative factor 
(DF), which has the high recognition performance when^ 
matching number between the input image and the right model 
is more than a matching number between the input image and 

25 the wrong model. 

Advantageous Effect 

The present invention has an effect to increase 
recognition performance of the iris recognition system and 
30 to reduce processing time for iris recognition, because the 
iris • recognition system can obtain an iris image 
appropriate for the iris recognition more effectively. 

The present invention detects a boundary between the 
pupil and the iris of an eye quickly and precisely, and 
35 extracts the unique features of the iris pattern by solving 

18 
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the problems of a non-contact iris recognition method, i.e., 
variation in the image size, tilting and moving, and 
detects a texture (iris pattern) by utilizing the Zernike 
moment having the visional recognition ability of a human 
5 being, regardless of motion, scale, illumination and 
rotation. . 

In the present, invention, an object in the iris image 
is identified by matching the feature quantity between 
models reflecting the reference value based on stability of 
10 the Zernike moment. and similarity between feature 
quantities in probability, and the iris recognition is 
verified by combining the LS and the Lmeds methods, to 
thereby authenticate the iris of the human being in rapid 
and precise. 

15 

Brief Description of the Drawings 

The above and other objects and features of the 
present invention will become apparent from the following 
20 description of the preferred embodiments given in 
conjunction with the accompanying drawings, in which: 

Fig. 1 is a block diagram showing an apparatus for 
extracting an iris feature and a system using the same in 
accordance with an embodiment of the present invention; 
25 Fig. 2 is a detail block diagram showing an apparatus 

for extracting an iris feature of Fig. 1 in accordance with 
an embodiment of the present invention; 

Fig. 3 is a flowchart describing, a method for 
extracting iris feature and a method for recognizing an 
30 iris using the same in accordance with an embodiment of the 
present invention; 

Fig. 4. is a diagram showing an appropriate iris image 
for the iris recognition; 

Fig. 5 is a diagram showing an inappropriate iris 
35 image for the iris recognition; 
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Fig., 6 is a flowchart showing a process for selecting 
an image at an image capturing unit in accordance with an 
embodiment of the present inventions- 
Fig. 7 is a graph showing a process for detecting an 
5 edge by using a 1-order differential operator in accordance 
with an embodiment of the present invention; 

Fig. 8 is a diagram showing a process for modulating 
connection ^ number for thinning" •in accordance with an 
embodiment of the present invention; 
10 Fig. 9 is a diagram showing a feature rate * of 

neighboring pixels for connecting a boundary in accordance 
with an embodiment of the present invention; 

Fig. 10 is a diagram showing a process for determining 
a center, of the pupil in accordance with an embodiment of 
15 the present invention; 

Fig. 11 is a diagram showing a process for determining 
a radius of the pupil in accordance with an embodiment of, 
the present invention; 

Fig. 12 is graphs showing a curvature graph and a 
20 model of an image in accordance with an embodiment of the 
present invention; 

Fig. 13 is a graph showing a process for transforming 
the image by using a linear interpolation in accordance 
with an embodiment of the present invention; 
25 Fig. 14 is a graph showing a linear interpolation in 

accordance with an embodiment of the present invention; 

Fig. 15 is a diagram showing a process for 
transforming a Cartesian coordinates system into a polar 
coordinates system in accordance with an embodiment of the 
30 present invention; 

Fig. 16 is a graph showing a Cartesian coordinates in 
accordance with an embodiment of the present invention; 

Fig. 17 is a graph showing a plane polar coordinates 
in accordance with an embodiment of the present invention; 
35 Fig. 18 is a graph showing a relation of zero-crossing 

20 
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points of ' first and second derivatives in accordance with 
an embodiment of the present invention,; 

Fig. 19 is a graph showing . a connection of zero- 
crossing points in accordance with an embodiment of the 
5 present invention; 

Fig. 20 is a diagram showing^ structures of a node and 
a graph of a two-dimensional histogram in accordance- with 
an embodiment oi- the present invention; 
* Fig-. 21 is a diagram showing a consideration when a 

10 transcendental probability is given in accordance ■ with an 
embodiment of the present invention; 

Fig. 22 is a diagram showing a sensitivity of a 
Zernike moment in accordance with an embodiment of the 
present invention; 
15 Fig'. 23 is a graph showing first and second ZMMs of an 

input image on a two dimensional plane in accordance with 
an embodiment of the present invention; 

Fig. 24 is a diagram showing method for matching local 
regions in accordance with an embodiment of the present 
20 invention; 

Fig. 25 is a diagram showing a False Rejection Rate 
(FRR) and a False Acceptance Rate (FAR) according to- a 
distribution curve in accordance with an embodiment of the 
present invention; 
25 Fig. 26 is a graph showing a distance distribution 

chart of an iris for an identical person in accordance with 
an embodiment of the present invention; 

Fig. 27 is a graph showing a distance distribution 
chart of an iris for another person in accordance with an 
30 embodiment of the present invention; 

Fig. 28 is a graph showing an authentic distribution 
and an imposer distribution in accordance with an 
embodiment of the present invention; and 

. Fig. 29 is a graph showing a decision of Equal Error 
35 Rate (EER) in accordance with an embodiment of the present 

21 
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invention. . . * 

Modes for Invention 

5 The above and other objects and features of the 

present invention will become apparent from the following 

description^ and thereby one of ordinary skill in the art 

can embody the principles of the present invention and 
invent various apparatuses within the concept and scope of 

10 the present invention. In addition, if further detailed 
description on the related prior arts is determined to blur 
the point of the present invention, the detail description 
shall be omitted. • Hereafter, preferred embodiments of the 
present invention will be described in detail with 

15 reference to the drawings. 

Fig. 1 is a block diagram showing an iris recognition 
system in accordance with an embodiment of the present 
invention. 

The iris recognition system includes basically an 
20 illumination (not shown) , a camera for capturing an image, 
e.g., desirably a digital camera (not shown), and can 
operates in a computer environment having such as a memory 
and a central processing unit (CPU) . 

The iris recognition system extracts features of an 
25 iris of a person by using an iris feature extracting 
apparatus having an iris image capturing unit 11, an image 
processing/dividing (fabricating) unit 12 and an iris 
pattern feature extractor* 13, and the iris feature is used 
for a verifying process of the person at an iris pattern 
30 registering unit 14 and an iris pattern recognition unit 16. 

At an initial time, a user must store feature data of 
its own iris in an iris database (DB) 15 and the iris 
. pattern registering unit . 14 registers the feature data. 
When verification is required later on, the user is 
35 required to identify him by capturing the iris using a 

22 
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digital camera, and then the iris pattern recognition unit 
16 verifies the user. 

When, the iris pattern recognition unit 16 verifies, 
the. captured iris feature is compared to the iris pattern 
5 of* the user stored in the iris DB 15. When the 
verification is • successful, the user can use the 
predetermined services. When the verification is . fa.iled^ 
the user is decided as an unregistered person or an illegal 
service user. 

io Detail structure of the iris feature extracting 

apparatus is as followings. As shown in Fig. 2, the iris 
extracting apparatus includes an image capturing unit 21, a 
reference point detector 22, an inner boundary detector 23, 
an outer boundary detector 24, an image * coordinates 

15 converter 25, an image analysis region defining unit 26, an 
image smoothing unit 27, an image normalizing unit 28, a 
shape descriptor extractor 29 and a reference value storing 
unit 30 and an image recognizing/verifying unit 31. 

The image capturing unit 21 digitalizes and quantizes 

20 an inputted image, and achieves an appropriate image for 
the iris recognition by detecting an eye blink and a 
location of a pupil and analyzing a distribution of 
vertical edge components. The reference point detector 22 
detects any reference point of the pupil from the acquired 

25 image and to thereby detect an. actual center point of the 
pupil. The inner boundary detector 23 detects an inner 
boundary wherein the pupil boundary on the iris. The outer 
boundary detector 24 detects an outer boundary wherein the 
iris borders on a sclera. The image coordinates, converter 

30. 25 converts a Cartesian coordinates system of a divided 
iris pattern image into a polar coordinates system and 
defines an origin of the coordinates as a center of a 
circular pupil boundary? . The image analysis region 
defining unit 26 classifies analysis regions of the iris 

35 image in order to uses the iris pattern defined based on 

23 
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clinical experiences of the iridology. The image smoothing 
unit 27 smoothes the image by filtering the analysis region 
of the iris image based on scale space in order to clearly 
distinguish a brightness distribution difference between 
5 neighboring pixels of the image. The image normalizing 
unit 28 normalizes a low-order moment with a mean size, 
v^herein the low-order moment is used for the smoothen image. 
The shape descriptor extractor 29' generates a Zernike 
moment based on the feature point' extracted from -the scale 

10 space and the scale illumination and extracts a rotation- 
invariant and noise-resistant by using Zernike moment shape 
descriptor. Also, the reference value storing unit 30 
(i.e., the iris pattern registering unit 14 and the iris DB 
15 of Fig. 1) stores a reference value as a template form 

15 by comparing stability of the Zernike moment to a 
similarity of Euclid distance, wherein the image pattern is 
projected into 25 spaces. . 

The image analysis region defining unit 26 is -n-bt an 
element included in the process of the iris recognition. 

20 The image analysis region defining unit 26 is included in 
the figure for the reference and shows that the feature 
point is extracted based on the iridology. The analysis 
region means the analysis region of the image appropriate 
for recognizing the iris that does not includes an eyelid, 

25 eyelashes or any predetermined part of' the iris to be 
intercepted by the mirror reflection from an illumination. 

Accordingly, the iris recognition system extracts the 
feature of the iris of the specific person by using the 
iris feature extracting apparatus 21 to 29, and recognizes 

30 the iris image i.e., identifies the specific person by 
matching the feature quantity between the reference value 
(the template) and a model reflecting stability and 
similarity of the Zernike moment of the iris image at the 
image recognizing/verifying unit 31 (i.e., the iris pattern 

35 recognition unit 16 of Fig. 1). 
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In particular*, the inner boundary detector 23 and the 
outer boundary detector .24 detect two reference points from 
a light source of the illumination, i.e., desirably 
infrared, from the eye image, determine a candidate pupil 
5 boundary point, determine a pupil location and a pupil size 
by obtaining a radius of a circle and a center point of a 
circle that are close to the candidate pupil . boundary based 
on the candidate center point, and to thereby detect the 
pupil region in real-time. In other word, the inner 

10 boundary detector 23 and the outer boundary detector 24 
detect two reference points by using a infrared 
illumination from the eye image acquired by the iris 
recognition system, determine candidate edge points between 
the iris and the pupil of the iris image where a line 

15 crossing the two reference points intersects, determine the 
candidate edge points where a perpendicular line crossing 
the center point between the two candidate edge points 
intersects, determine, the pupil location and the pupil 
size by obtaining the radius and the center point of the 

20 circle that are close to the candidate edge points based on 
the candidate center point where the perpendicular lines 
crossing the center point between the neighboring candidate 
edge points intersects, and to thereby detect the pupil 
region. 

25 The shape descriptor detector 29 detects the shape 

descriptor which is invariant to motion, scale, 
illuminatipn and rotation of the iris image. The Zernike 
moment is generated based on the feature extracted from the 
scale 'space and the scale illumination and the shape 

30 descriptor which is rotation-invariant and noise-resistant 
by using Zernike moment is extracted based on the Zernike 
moment. The indexed similar iris shape group database can 
be implemented based on the shape descriptor, and therefrom 
the indexed iris shape * group having the iris shape 

35 descriptor similar with that of the query image can be 
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15 



searched. - , e * 

The shape descriptor extractor 29 extracts the shape 
descriptor based on the linear shape descriptor extraction 
method. Thus, a skeleton is extracted from the iris image. 
A line list is obtained by connecting pixels based on the 
skeleton. The line list normalized is determined as the 
shape • descriptor. The iris shape database indexed by a 
dissimilar shape descriptor can be implemented by measuring 
dissimilarity of the indexed similar iris shape group based 
on .the linear shape descriptor and therefrom the iris image 
matched to the query image can be searched. 

Features of the each element 21 to 31 of the iris 
recognition system will be described in detail hereinafter 
in conjunction with Fig. 3. 

The iris image for the iris recognition irtust include 
the pupil, the iris furrow outside of the pupil and the 
entire colored part of the eye. Because the iris furrows 
are used for the iris recognition, the color information 
does not need. Therefore, a monochrome image is obtained. 

If the illumination is too strong, the illumination 
may stimulate the user's, eye, result unclear features of 
the iris furrows and can not prevent reflected rays to 
occur*. Under the consideration of the above conditions, 
the infrared LED is desirable. 
25 A digital camera, using a CCD or COMS chip that can 

achieve the image signal, display the image signal and 
capture the image. The image captured by the digital 
camera is preprpcessed. 

For simple description of the iris recognition phases, 
30 at first, the iris area included the eye must to be 
, captured. The resolution of the iris image is normally 
from 320x240 to 640x480. If there are a lot of noises in 
the image, an acceptable result can not be obtained even 
though preprocessing is excellently performed. Therefore, 
35 image capturing is important. It is important to maintain 

26 
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that conditions of' neighboring environment are unchangeable 
with the time. It is indispensable to determine a location 
of the illumination so that an interference of the iris by 
the reflected light due to the illumination must be 

5 minimized. 

Phases of, extracting the iris area and removing the 
noise from the image are called as' preprocessing. The 
preprocessing is required for extracting accurate iris 
features and includes a scheme for, detecting an edge 

10 between the pupil and the iris, dividing the iris area and 
converting the divided iris area into adaptable coordinates. 

The preprocessing includes detail processing phases 
that evaluate a quality of the achieved image, selects the 
image and makes the image to be utilized. A process to 

15 analyze the preprocessed features and convert the feature 
into a code having certain information is a feature 
extraction phase. The code is to be compared or to be 
studied. At first, the scheme for selecting the image is 
described, and then the scheme for dividing the iris will 

20 be described. 

The image capturing unit 21 achieves the image 
appropriate for the iris recognition by using 
digitalization, i.e., sampling and quantization, and 
suitability decision, i.e., eye blink detection, pupil 

25 location detection and vertical edge component distribution. 
The image capturing unit 21 performs to determine whether 
the image is appropriate for the iris recognition. The 
detail description will be described as follows. 

First of all, a method for selecting an image to be 

30 used efficiently through a simple suitability decision 
phase among a plurality of captured images from a fixed 
focus camera as a- method for achieving the image in the 
iris recognition system will be described. 

For achieving the image by the digital camera using 

35 the CCD or CMOS chip, a plurality of the images are 
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inputted and preprocessed in a determined time. A method 
for determining a moving image frame ranking through a 
real-time image suitability decision instead of recognizing 
all input images is used. 
5 According to the above processes, the processing time 

is decreased and the recognition performance is increased. 
For selecting an appropriate image, pixel distribution and 
edge component ratio are used. 

The digitalization at steps S301 and S302 of the 2- 

10 dimensional signals from the input image will be described. 

The image data is expressed' as an analog value of z- 
axis on the 2-dimensional space, i.e., x-y axis. For 
digitalizing the image, a space region is digitalized, and 
then a gray-level is .digitalized. The digitalization of 

15 the space region is called as a horizontal digitalization, 
and the digitalization of the gray-level is called as a 
vertical digitalization. 

The digitalization of the space region enlarges time 
axis sampling of a one-dimensional time series signal to a 

20 sampling of two-dimensional axis. In other words, the 
digitalization of the space region expresses the gray-level 
of discrete pixels. The digitalization of .the space region 
determines the resolution of the image. 

The quantization of the image-, i.e., the 

25 digitalization of the gray-level is a phase for limiting 
the gray-level into the determined steps. For example, if 
the nximber of the steps for the gray-level is limited to 
256, the gray-level can be expressed from 0 to 255. Thug, 
the gray-level is expressed in 8 bit binary number. 

30 . The image which is appropriate* for the iris 

recognition shows features of the iris pattern clearly and 
includes the entire iris area in the eye image. The 
accurate decision for the quality of the achieved imacje is 
important factor to affect the iris recognition system 

35 performance. Fig. 4 is an example image of a qualified 
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image for the iris recognition that the iris patt^ern is 
clear and there is no interference by the eyelid or the 
eyebrow. 

Meanwhile, if all input image are automatically 
5 provided to the iris recognition system, a recognition 
failure occurs due to an imperfect image and a low- 
qualified image. .There are four cases of failure eye image 
causing recognition fail. 

The first case is that there is an eye blink as shown 

10 in Fig. 5 (a), the second case is that a part of the iris 
area is truncated because a center of the pupil is out of 
the center of the image due to a user' s motion as shown in 
Fig. 5 (b) and the third case is that the iris area is 
interfered by the eyelash as shown in Fig. 5 (c) • An 

15 additional case is that there are many noises in the eye 
image (not shown) . Most of above cases fails to recognize 
the iris. Therefore, images of above cases are rejected by 
preprocessing, and to thereby improve efficiency of 
processing and a recognition rate. 

20 As mentioned above, the decision conditions for the 

qualified image can be provided with three functions as 
follows (See Fig. 6) at step S303. 

1) Decision condition function Fl: eye blink detection. 

25 2) Decision condition function F2: location of a pupil. 

3) . Decision condition function F3: vertical edge 
component distribution. 

The input image is subdivided into M x N blocks, whicH 
30 are utilized for functions of each step, and Table 1 as 
below shows an example of counting each block when the 
input image is subdivided into 3X3. 

Table 1 
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B3 


B4 


B5 


B6 


87 


B8 . 


• 



Because an eyelid area is brighter than a pupil area 
and .an iris area 'generally, it is determined" to the eye 
blink image when the image satisfies Eq. 1 as followings. 
Thus, it is determined that the image is unusable (i.e., 
the eye blink detection) . 



Max 



= ^M^ M,^Mean{B,) Eq. .1 
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15 



20 



The pupil is the region that has the lowest pixel 
values. The pupil region is located in as much as the 
center, the probability that the whole iris area is 
included in the image is increased (i.e., the pupil 
location detection) . Therefore, as Eq. 2, the image is 
subdivided into M X N blocks, the block of which pixel 
average is the lowest is detected, and the weight is added 
according to the location of the block. The weight of the 
pixel is smaller and smaller apart from the center of the 
image . 



ScoreiLoM(B%w) 

LoM(B) = Lociionq/Mm(B^,,yBj^^^^) 

F2 = w 



Eq. 2 



There are many vertical edge components at the pupil 
boundary and the iris boundary in the iris image (i.e., the 
25 edge component ratio investigation) . Based on a location 
of the pupil detected by Sobel edge detector as Eq. 1, the 
vertical edge components of the left and right region of 
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the image are investigated and the component comparisons 
are performed in order to investigate that whether an 
accurate boundary detection is possible or not and the 
change of the iris pattern pixel value is not large due to 
a shadow in the iris area extracting process which, is the 
next step of the image acquisition. 



20 



25 



30 



F3 = 



© = 



Eq. 



10 Wherein^ L is a left region of the pupil location, R 

is a right region of the pupil location, Ev is a. vertical 

component .and E^ is a horizontal component. 

The sum of each decision condition function value 

indicates utilization suitability of the image recognition 
15 process (Refer to Eq. 4), and is the base for counting 

frames of a moving picture achieved during a specific time 

(suitability investigation) • 



Eq. 4 



Wherein, T is a threshold, and intensity of the 
suitability is controlled according to the threshold. 

Meanwhile, the reference point detector 22 detects a 
real center point of the pupil after detecting a reference 
point of the pupil from the achieved image by Gaussian 
blurring, at step S304 including blurring, edge soften and 
noise reduction. Edge Enhancing Diffusion (EED) at step 
S305, image binalization at step S306. Thus, the noise is. 
removed by the EED method using a diffusion tensor, the 
iris image is diffused by Gaussian blurring, and a real 
center of the pupil is extracted by Magnified Greatest 
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Coefficient method. The diffusion is used for aecreasing 
bits/pixel of the image in the binalization process. Also, 
the EED method is used for decreasing edges. Detail part 
of the image is removed by Gaussian blurring that is a low 

5 frequency pass filter. When the image is- diffused, the 
actual center and size of the pupil are found by changing a , 
threshold used for the binalization process. Detail' 
description is as follows. 

As the first preprocessing step, the edge is softened 

10 and noises in the image are removed by Gaussian blurring at 
step S304. However, too large Gaussian value cannot be 
used because dislocation occurs in the low-resolution image. 
If too large Gaussian deviation value is used, dislocation 
occurs in a low-resolution image. If there is mere noise 

15 in the image, Gaussian deviation value can be small or none. 

Meanwhile, at step S305, the EED method is applied 
strongly to a part where the direction is the same with the 
edge, and is applied weakly to a part where the direction 
is an orthogonal to the edge by considering the local edge 

20 direction. Non-linear Anisotropic Diffusion Filtering 
(NADF) is one of the diffusion filtering methods and the 
EED method is a major methocj of the NADF. 

In the EED method, the iris image after Gaussian 
blurring is diffused and a diffusion tensor matrix is used 

25 by considering not only a contrast of the image but the 
edge direction. 

At the first phase for implementing the. EED method, 
the diffusion tensor ' instead of a conventional scalar 
diffusivity i%.used. 

30 . The diffusion tensor matrix can.be calculated based on 

eigenvectors vl and v2. The vl is paralleled with Vu as Eq. 
5 and the v2 is orthogonal to Vu as Eq. 6. 

V m Vu . ^ ' Eq. 5 
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^2 VU Eq^ 6 

Therefore^ Eigenvalues Kl and A.2 are selected in; order 
5 to perform smoothing at the part paralleled with the edge 
rather than the part orthogonal to the edge. Eigenvalues 
are expressed as: 

diffusion across edge \^ '= D( I Vu \^2) ^ 
diffusion along edge X;2:= 1 g 



10 



15 



According to the above method, the diffusion tensor 
matrix D is calculated based on an equation expressed as: 

In order to implement the diffusion tensor matrix D in 
a real program, the vl and v2 must be clearly defined. If 

20 the original iris image is expressed as Gaussian-filtered 
vector (gx, gy) , the vl makes the original iris image to be 
a parallel with Gaussian filtered image and can be 
expressed as (gx, gy) as shown in Eq, 5. The v2 is 
orthogonal to Gaussian-filtered image and the scalar 

25 product of (gx, gy) and the v2 is made to be zero as shown 
in Eq. 6. Therefore, the v2 is expressed as (-gx, gy) . 
Because vl' and v2' are transpose matrix of the vl and the 
v2 respectively, the diffusion tensor matrix* D can be 
expressed as: 
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Eq. 10 



Wherein the d can be calculated based on dif f usivity 
which is presented in Eq. 14 as follows. 
5 At the second phase for implementing the EED method^ a 

constant K is determined. The K denotes how much an 
absolute value is accumulated in a histogram of the 
absolute value. If the K is 90* or above, it can be a 
problem that detail structures of the iris image is quickly 

10 removed. If the K is 100% , it can be a problem that the 
whole iris image is blurred and the dislocation occurs. If 
the K is too small, the detail structures still remain 
after a lot of time iterations. 

At the third . phase for implementing the EED method, 

15 the diffusivity is evaluated. A gradient is calculated by 
Gaussian blurring the original iris image. A magnitude of 
the gradient is obtained. Because a gray-level is rapidly 
changed at the edge, a differential operation that takes 
the gradient is used for extracting the edge. The gradient 

20 at point (x, y) of the iris image f(x, y) is a vector 
expressed as Eq. 11. The gradient vector at point (x, y) 
denotes maximal change rate direction of the f. 




25 

The gradient vector Vf is expressed as: 

r 
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The Vf is equal to a maximal increase rate per unit 
length at a direction of Vf. 

In practice, the gradient is approximated as shown in 
5 Eq. 13 expressed with absolute values of the gradient* Eq. 
13 is easy to calculate and implement with a limited 
hardware . 



10 



15 



Eq. 13 



The diffusivity expressed as Eq. 14 is obtained based 
on the K and the obtained at the second phase. 

D = 1 / (1 + magnitude of gradient/ K^2) 



At the forth phase for implementing the EED method, 
the diffusion tensor matrix D is obtained as shown in Eq. 
10 and a diffusion equation is "evaluated based on Eq. 15. 
At first, the gradient of the original iris image and then 
20 the gradient of the Gaussian-filtered iris image are 
applied to the original iris image. For the gradient of 
the Gaussian-filtered iris image does not exceed 1^ the 
normalization must be performed. 

25 at u =d iv (D-Vu) Eq. is 

The iris image is dif.fused under consideration of not 
the edge direction but . contrast because the diffusion 
tensor matrix is used. The smoothing is weakly performed 
30 where orthogonal to the edge, and is strongly performed 
where paralleled with the edge. Therefore, a problem that 
the edge with the noise is extracted where there are a lot 
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of noises in the edge can be improved. i 

A process from the second phase to the forth phase is 
j-epeated up to the maximal time iteration. Problems caused 
by many noiises in the original iris image, scale-invariant 
5 image due to the constant K and unclear edge extraction due 
to noises at the edge are solved by processing the aboye 
four phases. 

The Vu as shown in Eqs. 5 to 15 denote the diffusion 
of each part of the image*. The diffusion tensor matrix D 
10 is evaluated based on the eigenvector for the edge of the 
image and then the divergence is performed resulting linear 
integral, and thereby the contour of the image is obtained. 

Meanwhile, the iris image is transformed into a binary 
image for obtaining a shape region of the iris image at 
15 step S306 (the image binalization) . The binary image is 
black and white data of the monochrome iris image based on 
the threshold value. 

For image subdivision, gray-level or chromaticity of 
the iris image is evaluated into the threshold value. For 
20 example, the iris area is darker than the retina area of 
the iris image. 

Iterative thresholding is used for obtaining the 
threshold value when the image binalization is performed. 

The iterative thresholding method is to improve an 
25 estimated threshold value by the iteration. It is .assumed 
that the binary image . obtained based on the first threshold 
•is used for selecting the threshold resulting a better 
image. A process for changing the threshold value is very 
important to the iterative thresholding method. 
30 At the first phase of the iterative thresholding 

method, an initial estimated threshold value T is 
determined. A mean brightness of the binary image can be a 
good threshold value. 

At . the second phase of the iterative thresholding 
35 method, the binary image is subdivided into a first region 
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10 



Rl and a second region R2 based on the .initial estimated 
threshold value t. 

At the third phase of the iterative thresholding 
method, average gray levels yl and \i2 of the first region 
Rl and the second region R2 are obtained. 

At the forth phase of the iterative thresholding 
method, a new threshold value is determined based on Eq, 16 
expressed as: 

Eq. 16 



At the fifth phase r a process from the second phase to 
the forth phase is iterated until the average gray levels 
pi and \i2 are not changed. 

After the binalization of the whole image, data is 
15 obtained. The inner boundary and the outer boundary are 
detected based on the data. 

A process for detecting the inner boundary and the 
outer boundary is described as follows, i.e., a pupil 
detection that determines a center and a radius of the edge 
20 at steps S307 to S309. 

The inner boundary detector 23 detects the inner 
boundary between the pupil and the iris at steps S307 and 
S308. The binary image binalized based on Robinson compass 
Mask is subdivided into the iris and the background, i.e., 
25 the pupil. And, intensity of the contour is detected based 
. on Difference of Gaussian (DoG) so that only intensity of 
contour is appeared. Then, thinning is performed on. the 
contour of* .the binary image using Zhang Suen algorithm. 
The center coordinate is obtained based on bisection 
30 algorithm. A distance from the center coordinate to a 
radius of the pupil in the counter clock-wise is obtained 
based on Magnified Greatest Coefficient method. 

The Robinson compass Mask is useci for detecting the 
contour. The Robinson compass Mask is a first-order 
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differentiation and a form of 3 X 3 matrix that evaluates 
an 8-<iirectional edge mask by rotating Sobel edge sensitive 
' a diagonal directed contour to the left. 

Also^ the DoG that is a quadratic differentiation is 
5 used for extracting the detected contour. The DoG 
decreases noises in the image based on Gaussian smoothing 
- fufiction, decreases a lot of the operations due to a mask 
size by decreasing two Gaussian mask, i.e.> LoG, and is a 
high frequency pass filtering operation. The high 

10 frequency denotes that a brightness distribution difference 
with the background is large. Based on the above 
operations, the contour is detected. 

Also, the thinning transforms the contour into a line 
of one pixel and obtains the center coordinate using the 
15 bisection algorithm, and to thereby obtain the radius of 
the pupil based on Magnified Greatest Coefficient method. 
The contour is formed to a circle and then^ the center 
point is applied to the circle^ and thereby the most 
similar shape to the pupil is selected. 
20 The outer boundary detector 24 detects * the outer 

boundary between the iris and the sclera at steps S307 to 
S309. 

For the outer boundary detection^ the center point is 
obtained based on the bisection algorithm. A distance from 

25 the center point to a radius of the pupil is obtained based 
on Magnified Greatest Coefficient ^method. Wherein, the 
linear interpolation is used to prevent that the image is 
distorted when coordinates system is ' transformed from 
(^artesian coordinates system to the polar coordinates 

30 system. " - ^ 

Edge extraction of the image, i.e., thinning and 
labeling, is needed at step S307 for the inner boundary and 
the outer boundary detections at steps S308 and S309. The 
edge extraction of the image means a ^process that the 

35 binary image is subdivided into the iris and the background 
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based on- the Robinson compass Mask, the intensity of the 
contour is enhanced based on the DoG, and the thinning is 
performed on the contour based on the Zhang Suen algorithm. 
.Referring to the edge extraction at step S307, because 
5 the edge is where the density is rapidly changed, the 
differentiation analyzing the value of the function change 
is used to extracting the contour. There . 'are a* first 
differentiation, i.e., the gradient and a quadric 
differentiation, i.e., the laplacian in the differentiation. 
10 Also, there is the edge extracting method by using a 
template-matching. 

The gradient observes a brightness change of the iris 
and is a vector G(x, y) = (fx, fy) expressed as: 



15 



20 



25 



Wherein, the fx is a gradient of x direction and the 
fy is a gradient of y direction. 

The Robinson compass Mask gradient operator 3X3 is 
illustrated in below and is the 8-directional edge mask 
made by rotating the Sobel mask to the left. The direction 
and the magnitude are determined according to the direction 
and the magnitude of the mask having the maximum edge value. 



The contour of the image must be pre-extracted to 
preprocess the acquired image. The iris and the background 
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are , subdivided based on the Robinson compass Mask that is 
the gradient- The gradient at the point (x, y) of the 
image f(x^ y) is expressed as Eq. 18. A magnitude "of the 
gradient vector (Vf) is expressed as Eq. 19. The gradient 
5 based on the Robinson compass Mask is given from the 
maximum edge mask among the following 8-directional masks 
based on Eq. 20. The z is brightness of pixel overlapped' 
by the mask at a location. The edge direction is a 
direction where the edge is put and can be cierived from a- 

10 result of the gradient. The edge direction is orthogonal 
to the gradient direction. That is/ the gradient direction 
denotes by a direction where difference value is changed 
largely and the edge must exist where the valued is changed 
largely. Therefore the edge is orthogonal to the gradient 

15 direction. 

Fig. 7 (b) is an image having the extracted contour. 



\7F 



dx 



Eq. 18 



20 ^ I «l I ^, 

<7, = (Z, + Z» + Z^) - <Z, + 2Za + Z,) 

G, =(Z, + Z^ + Zs,)-(Z, +2Z^ +Z,) Eq. 20 

The subscripts denote pixels as shown in Eq. 20. 
25 Meanwhile, the laplacian observes the brightness 

distribution difference with neighboring area. The 
laplacian performs the differentiation on the result of the 
gradient, and to thereby detect the intensity of the 
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contour. That is, only magnitude of the edge but not tHe 
. direction is obtained. The lapla'cian operator targets to 
find zero-crossings where the value is changed from + to - 
or from - to +. The laplacian decreases the noise in the 
5 image based on the Gaussian smoothing function and uses the 
DoG operator mask that decreases many operations due to the 
mask magnitude by subtracting the Gaussian masks having 
different values. Because the DoG approximates the LoG, a 
desirable 'approximation is obtaine^d. when a ratio 01/02 is 
10 1.6. 

The LoG and the DoG of two-dimensional function f(x, 
y) are expressed as: 



Eq. 21 



15 



Eq. 22 



The edge detection using the laplacian operator uses 
the 8-directional laplacian mask as shown in Eq. 23 and 8 
20 direction values based on the center, and to thereby 
determine a current pixel value. 

• ' 

; Loplac^»n(x.y> « gxfCx.y) - { K^c .y- 1) Kx . y+ |> + f(x-l. y> ^ f<x* I .y)* r6c* 1 .y^ 1) 
f(x-l,y-0-»-.f<jc-l,y-^ fCx^l.y-1)} - 

Eq. 23 

25 

The laplacian quadric differentiation operator 3X3 is 
as followings. 

Laplacian mask: direction-invariant 
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The thinning is described hereinafter. 

The Zhang Suen thinning algorithm is one of parallel 
5 processing" methods, wherein deletion means that a" pixel is 
deleted for , the thinning. Therefore, the black is 
converted into the white. 

Connection number is a nxamber indicating whether a 
pixel is connected to neighboring pixels or not. That is, 

10 if the connection number is 1, the center pixel, i.e., 0, 
can be deleted. A convergence from black to white or from 
white to black is monitored. Fig. 8 shows a check all 
pixels are converted from the back to the white. The pixel 
must be 1 regardless neighboring pixel numbers. 

15 Meanwhile, a labeling means distinguishing iris 

sessions apart from each other. A set of neighboring 
pixels is called as a connected component in a pixel array. 
One of most frequently used operations in a computer vision 
is to search the connected component from the given image. 

20 Pixels belongs to the connected component have high 
probability to indicate an object. A process for giving 
the label, i.e., the number, \o the pixels according to the 
connected component where the pixels belongs is called as 
the labeling. An algorithm for searching all connected 

25 components, giving the same label to pixels included in an 
identical connected component is called as a component 
labeling algorithm. The sequential algorithm takes short 
time and small memory comparing to an iteration (algorithm, 
and completes calculations within two times scanning to the 

30 given image. 

The labeling can be completed with two loops using an 
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equivalent table. . The drawback is that the labeling 
numbers are not continuous. The entire iris sessions .are 
checked and labeled. During the labeling, if other label 
. is detected^ the label is inputted in the equivalent table. 
5 The labeling is performed with the minimum label in a new 
loop. 

At first, a black pixel on the boundary is searched as 
shown in .Fig. 9. The boundary, point has 1-7 white pixels 
in the neighbor based on a center pixel. An isolate point 

10 is excluded. The isolated point's all neighboring pixels 
are black. Then, the labeling is performed in a horizontal 
direction and then a vertical direction. With two 
directional labeling as above, a U shape curve can be 
labeled in one time, and thereby the time is saved. 

15 The center point of the boundary and the radius 

determination, i.e., the pupil detection, steps for the 
pupil detection at the inner boundary detector 23 and the 
outer boundary detector 24 will be described. 

As above mentioned, in the pupil detection process, 

20 two reference points of the pupil from the light source of 
the infrared illumination are detected at SI. The 
candidate boundary points are determined at S2. The pupil 
region is detected in real-time by obtaining the radius and 
the center point which are the closet to the candidate 

25 boundary point based on the candidate center point and 
determining the pupil location and the pupil size at S3. ^. 

The process for detecting two reference points in the 
^ pupil from the light source of. the infrared illumination 
will be described. ^ _ _ . . 

30 For detecting the pupil location, the present 

invention obtains a geometrical variation of the light 
component generated in the eye image, calculates an average 
of the geometrical variation and uses the average as a 
template by modeling the average into the Gaussian waveform 

35 as Eq. 24. 
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. ■ V. ■ ' , Eq. 24 

Wherein, x is a horizontal location, y is a vertical 
5 location and O is a filter size. 

The two reference points are detected by performing a 
template matching based on the template so that, the 
reference point is selected in the pupil of the eye image. 

Because the illumination in the pupil of the eye image 
10 is the only part where a radical change of the gray-level 
occurs, it is possible to extract the reference point 
stably. 

The process for determining the • candidate pupil 
boundary point at S2 is described hereinafter. 

15 At the first step, a profile is extracted presenting 

the pixel value change of the waveform in +/- x axes based 
on the two reference points. The candidate boundary masks 
h(l) and h(2) corresponding to the gradient are generated 
in order to detect two candidate boundaries passing the two 

20 reference points in form of one-dimensional signal in the x 
direction. Then, the candidate boundary point is 

determined by generating a candidate boundary waveform (Xn) 
using a convolution of the profile and the candidate 
boundary mask. 

25 At the second step, another candidate boundary • point 

is determined by the same method of the first step on a 
perpendicular line based on the center, point, bisecting a 
distance between the two candidate boundary points. 

Meanwhile, the process for detecting the pupil region 

30 in real time by obtaining the radius and the center point 
of a circle which are the closet to the candidate boundary 
point based on the candidate center point and determining 
the pupil location and the pupil size at S3 will be 
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described hereinafter. 

The radius and the center point of the. circle closet 
to the candidate boundary point is obtained by using the 
candidate center point where the perpendicular lines at the 
5 bisecting points between the neighboring candidate boundary 
points are intersected. Hough transform for obtaining a 
circle component shape is applied to the above method. 

It is assume that there are two points A and B on a 
circle and a point C is "a bisecting point of a* line AB 
10 connecting points A and B. A line that crosses the point C 
and is perpendicular to the line AB always passes an origin 
O of the circle. An equation of a line OC is expressed as: 

j; = - X + ^/^ ^^'v" Eq. 25 

15 

In order to obtain the features and the location of 
the connected components group that make the circle, the 
center point is used as an attribute of the connected 
components group. Because the center of the inner boundary 

20 of the iris is changed and the boundary is interfered by 
the noise, a conventional method for obtaining a circle 
projection may evaluate an inaccurate pupil center. 
However, because the method uses the two light sources that 
are apart from a specific distance, the candidate center 

25 distribution coefficient of the ^bisecting perpendicular 
lines i§ appropriate to determine the center of the circle. 
Therefore, a point where the perpiendicular lines are mostly 
crossed among the candidate Isenter points is determined as 
.the- center 'of the circle (See Fig. iV. ) . 

30 After extracting the center of the circle according to 

the above method, the radius of the pupil is determined. 
One of the radius decision methods is an average method. 
The average method is to obtain an average distance of all 
distance of the group components making the circle from the 
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determined center point. That is similar' with Daugmans' 
method and Groen's method, if there are many noises in the 
-.image, the circumference component is distortedly 
recognized and the "clistortipn affects to the pupil radius. 
5 With comparison to the above method. Magnified 

Greatest Coefficient method is based oh the enlargement 
from a small region to a large region. At the first step, 
a longer distance is selected among pixel distances between 
the center point and the candidate boundary points. At the 

10 second step, the range becomes narrower by applying the 
first step at the candidate boundary points over the 
selected distance. therefore, the radius representing the 
circle is obtained by searching an integer finally. 
Because the distribution of transformation in all 

15 directions due to a contraction, expansion and a horizontal 
rotation of an iris muscle must be considered when the 
above method is used, it can extract the inner boundary of 
a stable and identical iris region (See Fig. 11.) 



20 Eq. 26 

Coordinates of the y is determined based on the radius 
and Eq. 26. If there is the black pixel in the image, the 
center point is accumulated. The circle is found based on 

25 the center point and the radius by searching the maximum 
accumulated center point. (Magnified Greatest Coefficient 
method) - 

The center point is obtained using a bisection 
algorithm. Because the pupil has different curvature 

30 according to the kind, the radius is obtained based on the 
Magnified Greatest Coefficient method in order to measure 
the curvature of the pupil. Then, a distance from the 
center point to the outline in the ' counter clock-wise is 
obtained. It is presented on a graph that an x-axis is a 
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rotation angle and a y-axis is a distance from- the center 
to the contour. In order to find the features of the" image, 
a peak and a valley of the curvature are obtained and the 
maximum length and an average length between the curvatures 
5 are evaluated. 

Fig. .12 is a graph showing the curvature graph of the 
acquired circle image (a) and the acquired star-shaped 
image (b) . In the case of the circle image (a) , because 
the distance from the center to the contour is uniform, the 

10 y has fixed value and the peak and the valley are r. The 
above case is weak for drape property. If the image is 
drifted, the distance from the center to the contour is 
changed. Therefore, the y is changed and has the curvature 
in the graph. In the case of the star-shaped image (b) , 

15 there are four curvatures in the graph, and the peak 
becomes r and the valley becomes a. 

Circularity shows how much the image looks like a 
circle. If the circularity is close to 1, the drape 
property is weak. If the circularity is close to 0, the 

20 drape property is strong. For evaluating the circularity, 
a circumference and an area of the image are needed. The 
circumference of the image is a sum of distances between 
pixels on the outer boundary of the image. If the pixel of 
the outer boundary is connected perpendicularly or in 

25 parallel, the distance between pixels is 1 unit. If the 
pixel is connected in diagonal, the distance between pixels 
-is 1.414 units. The area of the image is measured as a 
' total number of the pixels inside of the outer -boundary. A 
formula for obtaining the circularity is expressed as: 

30 - -. 

circularityie) = ^ — ^^rycarea ^ 27 

[circumference) 

According to the edge extraction process at step S307, 

the inner boundary is confirmed, and the actual pupil 
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center is obtained using the bisection algorithm. ,Then, 
the radius is obtained , using the Magnified. Greatest 
Coefficient method wKen the pupil is assumed to be a 
perfect circle/, and the distance from the center to the 
5 inner boundary in the counter clock-wise is measured^ and 
thereby the data is generated as shown in Fig* 12 (the 
inner boundary detector 23 and the outer boundary detector 
24 perform. ) 

The processes from the binalization at step S306 to 
10 the inner boundary extraction at step S308 are summarized 
in sequence as follows: EED -> binalization of -> edge 
extraction -> bisection algorithm -> Magnified Greatest 
Coefficient method -> inner boundary data generation -> 
image coordinates system transformation, 
15 Meanwhile, in the outer boundary detection at step 

S309, the edge betweien the pupil and the iris is found with 
the same method of the inner boundary detection filtering, 
i.e., the Robinson compass mask, the DoG and the Zhang Suen. 
Wherein, where the difference between the pixels is a 
20 maximum is determined as the outer boundary. The linear 
interpolation is used in order to prevent that the image is 
distorted due to motion, rotation, enlarge and reduction 
and in order to make the outer boundary as the circle after 
thinning. 

25 The bisection algorithm and the Magnified Greatest 

Coefficient algorithm are used in the outer boundary 
detection at step S3G9. Because the gray-level difference 
of the outer boundary is riot clearer than that of the inner 
boundary, the linear interpolation is used. 

30 The process of the outer boundary detection at step 

S309 is described hereinafter. 

Because the iris boundary is blurred and thick, it is 
hard to find the boundary exactly. The edge detector 
defines where the brightness is changed most as the iris 

35 boundary. The center of the iris can be searched based on 
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the pupil center, and the iris radius can be searched based 
on that the ir>s radius is mostly uniformed .in the fixed 
focus came;ra • 

The ^dqe between the pupil, and the iris . is obtained 
5 with the same method of the inner boundary detection 
filtering, and where the pixel difference is a maximum is 
detected as the outer boundary by checking the' pixel 
difference • 

Wherein, the transformation, i.e., motion, rotation, 
enlargement and reduction, using the liner interpolation is 
used (See Fig. 13.) 

As shown in Fig. 13, because a pixel coordinates is 
not matched 1 to 1 if the image is transformed,- the inverse 
transformation complements the above problem. Wherein, if 
there is a pixel that is not matched in the image, the 
pixel is shown based on the pixel of the original image. 

The linear interpolation as shown in Fig. 14 
determines a pixel based on four pixels based on how close 
X, y coordinates are. 
20 It is expressed using p and q as: 

p(q* equation+ ( 1 - q) * eqiiation) + q(p* equation+ (1 - p) * equation) . 

The image distortion is prevented by using the linear 
interpolation . 

25 The transformation is subdivided into three cases, 

i.e.^ motion, enlargement & reduction and rotation. 

The motion is easy to transform. A regular motion is 
to subtract a constant and n inverse motion is to plus the 
constant expressed as: . 



30 



Eq. 28 



The enlargement is to divide by the constant as Eq. 29 
below. Therefore, x and y are ' enlarged. Also, the 
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reduction is to multiply the constant, 



x/ a 



Eq.- 29 



5 The rotation is to use a rotation transformation 

having a sine function and a cosine function expressed as: 



SirtO CosO O 
Cc?sO — SinO O 
O O 1 



Eq. 30 



10 By unfolding Eq. 30, the inverse transformation 

equations are derived expressed as: 

15 The processes from the binalization at step S306 to 

the outer boundary extraction at step S309 are summarized 
as follows: EED -> iris inner /outer binalization -> edge 
extraction -> bisection algorithm -> Magnified Greatest 
Coefficient method -> iris center search -> iris radius 

20 search -> outer boundary data generation -> image 
coordinates system transformation. 

The process for transforming the Cartesian coordinates 
system into the polar coordinates system at step S310 will 
be described. As shown in Fig. 15, the divided iris 

25 . ^pattern image is transformed from the Cartesian coordinates: 
system into the polar coordinates system. The divided iris 
pattern means a donut-shaped iris. 

The iris muscle and the iris layers reflect a defect 
of the structure and the connection state. Because the 

30 structure affects to a function and reflects the integrity. 



50 



wo 2005/024708 PCT/KR2004/00228S 

the structure indicates the resistance of the organic and 
the genetic stam|>. The related signs are Lacunae, Crypts, 
Defect signs and Rarifition. 

.In order to use the iris pattern based on the clinical 
5 experience of the iridplogy as the features, the image 
analysis region defining unit 26 divides the iris analysis 
region as follows. Thus, it is subdivided into 13 sectors 
based on the clinical experience of the iridology. 

Therefore, the region is subdivided into a sector 1 at 

10 right and left 6 degree based on the 12 clock direction, a 
sector 2 at 24 degrees, in the clock-wise, a sector 3 at 42 
degree, a sector 4 at 9 degree, a sector 5 at 30 degree, a 
sector 6 at 42 degree, a sector 7 at 27 degree, a sector 8 
at 36 degree, a sector 9 at 18 degree, a sector 10 at 39 

15 degree, a sector 11 at 27 degree, a sector 12 at 24 degree 
and a sector 13 at 36 degree. Then, the 13 sectors are 
subdivided into 4 circular regions based on the iris. 
Therefore, each circular region is called as a sector 1-4, 
a sector 1-3, a sector 1-2, a sector 1-1, and so on. 

20 Wherein, 1 sector means Ibyte and stores iris region 

comparison data in a parted region, and to thereby be used 
for determining the similarity and the stability. 

The two-dimensional coordinates system is described as 
follows . 

25 The Cartesian coordinates system is a typical 

coordinates system showing 1 point on a plane as shown in 
Fig. 16. A point O is determined as origin on the plane, 
and two perpendicular lines XX'' and YY' crossing origin are 
axes. A point P on the plane is presented with a segment 

30 OP' = X passing the point P and parallel with' x-axis and 
with a segment OP' ' = y passes passing the point P and 
parallel with y-axis. Therefore, the location of the point 
P is matched to an ordered pair of two real numbers (x, y) , 
and reversely the location of the point P can be determined 

35 from the ordered pair (x,y). 
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^ Plane polar coordina1:es system ''is a coordinates system 
presented with a length of a segment connecting a point on 
the plane and the origin and an angle of the segment and an 
axis passing the origin. A 'polar angle 9 has a plus, value 
5. in the counter clock-wise of the mathematical coordinates 
system^ but the polar angle 9 has a plus value in the 
clock-wise * of the general measurement such as the azimuth 
angle. 

Referring to Fig. 17, the 9 is a polar angle, the O is 
10 -a pole, and the OX is a polar axis. 

A relation of the Cartesian coordinates system (x, y) 
and the Plane Polar Coordinates system (r, 9) is expressed 
as: 

r=rx* -H y^, e=tan~^(y/x) 

X = r cos e. y =r sine 3^ 

The image smoothing at step S311 and the image 
normalization at step S312 will be described. 

The image normalizing unit 28 normalizes the image by 
20 a mean size based on a low-order moment at step S312. 
Before the normalization, the image smoothing unit 27 
performs the smoothing on the image by using Scale-space 
filtering at step S311. When a gray-level distribution of 
the image is weak, it is improved by performing a histogram 
25 smoothing. Therefore, the .image smoothing is used for 
distinguishing clearly the gray-level distribution, 
difference among neighboring pixels. The scale-space 
- - filtering is performed in the .image smoothing process. The 
Scale-space filtering is a form that Gaussian function and 
30 the scale constant is combined, and is used for making a 
size-invariant Zernike moment after the normalization. 

The normalization at step 3312 and then the image 
smoothing at step 8311 will be described. 
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The . normalization at. step S312 . must be performed 
before a post processing is performed. The normalization 
uniforms . the size of the image, defines locations and 
adjusts a thickness of the line, ^nd to thereby standardize 
5 the iris image. 

Th^ iris image can be characterized based on 
topological features. The topological feature is defined 
as invariant features in spite of (elastic deformations of 
the image. Topological invariance excludes connecting 

10 other regions or dividing other regions. For a binary 
region. Topological characteristic features include the 
number of the hole and embayment, protrusion. 

More precise expression than the hole is a subregion 
which exists inside of the iris analysis region. The 

15 subregion can appear recursively. The iris analysis region 
can include the subregion including another subregion. A 
simple example for explaining a discrimination ability of 
Topology is an alphanumeric. Symbols 0 and 4 have one 
subregion and B and 8 have two subregions. 

20 Evaluation of the moments indicates a systemic method 

of the image analysis. The most frequently used iris 
features are calculated based on three moments from the 
lowest order. Therefore, the area is given by 0-order 
moment and indicates the total number of the region-inside. 

25 A centroid determined based on 1-order moments provides the 
measurement of the shape location. A directional motion of 
the regions is determined Ipased on principal axes 
determined by order moments. 

Information- of the low-order moments allows evaluating 

30 7central moments, normalized central moments and moment 
invariants. These quantities delivery shape features that 
are invariant to the location, the size, and the rotation. 
Therefore, when the location, the size and the directional 
motion does not affect to the shape identity, it is useful 

35 for shape recognition and the matching. 

53 



wo 2005/024708 



PCT/KR2004/002285 



•The moment analysis' is based on the pixels inside of 
the iris shape region* Therefore^ a growing or a filling 
of the iris shape region for summing " all pixels inside of 
the iris shape region are needed in advance. The moment 
5 analysis is based on the contour of the bounding region of 
the iris shape image, and it requires the. contour detection. 

The pixels of the bounding region are allocated as 1 
(ON) for the binary image of the iris analysis region, and 
■ a moment mpq of the binary, image is defined as Eq. 33 below: 

10 Region-Based Moments, (p+q) ^^-order normalized moments for 
2-order iris analysis region shape f (x, y) are expressed as 
Eq. 33 below. Wherein, when p=0 and q=0, 0-order 
normalized moment is defined as Eq. 34 below and indicates 
a pixel number included in the iris analysis region shape. 

15 Therefore, the measurement of the area is provided. 
Generally, the number of the shape indicates a size of the 
shape but is affected by the threshold value in the 
binalization. Even though the size of the shape is same, 
the contour of the iris image resulted by the binalization 

20 based on a low value is thick and the contour of the iris 
image resulted by the binalization based on a high value is 
thin. Therefore, the pixel number varies in large at the 
0-order moment value. 



25 



«=0^ Eq. 33 



N M 



«=0.v=0 Eq. 34 



Table 2 
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1^- . . ■ 



' 1^11 

1^ 



[Moments and Vertex Coordinates] 

Generally/ the moitient is defined based on the pixel 
5'. • location- and the pixel value expressed" as: 

ntp^ = £^x^y^f{x,y)dx(fy Eq. 35 



Moment equations up to the quadratic-order are easily 

10 derived based on a static point defining the bounding 
region contour of the binary iris shape simply connected. 
Therefore, if it is possible to express a polygonal of a 
region contour, the area centroid and the directional 
motion of the principal axes can be easily derived based on 

15 the equation in Table. 2. 

* the lowest-order moment moo indicates the total pixel 
number inside of the iris analysis region shape and 
provides the measurement of the area. If the iris shape in 
the iris analysis region is particularly larger or smaller 

20 than another shape in the iris* image, the lowest -order 
moment m^o is useful as the shape descriptor. However, 
because the area occupies smaller part or larger part of 
the shape according to the scale of the image, a distance 
between the object and the observer and a perspective, it 

25 can not be used imprudently. 
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The 1-order moment of the x and the y normalized based 
on the area of the iris image provides coordinates of the x 
^fiX and. y centroid. The average location of the iris shape 
rigion is determined based on the coordinates of the x and 
5 y centroid. 

After the * iris shape division process, all shapes of 
the image are given the same label and then the up and down 
boundaries of the iris are denoted by A and B, the left and 
right boundaries of the iris are denoted by L and R 
10 respectively, and the coordinates of the x and y centroid 
are expressed as: 

Eq. 36 



15 



The central moment indicates iris shape region 

descriptor normalized based on the location. 



/'«=Z(^-*cyO'->'c)' Eq. 37 

R 

20 

Generally, central moment is normalized with 0-order 
moment as Eq. 38 in order to evaluate the normalized 
central moment. 

25 Vpg=Mpg/Mm ,r = (j^ + ^)/2 + l Eq. '38 

_ _ ._ _ ^ 

The; normalized central moment which is the most 
frequently used is a M n that is a 1-order central moment 
between x and y. The |J ^ provides the measurement of the 
30 variation from the circle regions shape. Therefore, a 
value close to 0 describes a region similar to the circle 
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and a large value describes a region dissimilar to the 
circle, A principal major axis is defined as an axis 
passing the centroid having the maximum inertia moment and 
a principal minor axis is defined as an axis passing the 
minimum centroid. Directions of the principal major and 
minor axes are given as: 



tan^ = - 
2 



Eq. 39 



10 Estimation of the direction provides an independent 

method for determining an orientation of an almost circle 
shape. TherefpfeV. it . is. an appropriate .par:aineter. to. 
monitor the orientation motion . of the. transformed corltour/ 
e.g., for time-variant shapes. 

15 The normalized and central normalized moments . are 

normalized based, on the scale (area) and the motion 
(location) . The normalization based on the orientation is 
provided by a family of the moment invariants. Table 3 
evaluated based on the normalized central moments shows 

20 four moment invariants from the first. 



25 



Table 3 

Central moments 



30 



Moment invariants 

^1 =720 +7o2 y ^2 =(720 -702)^ -^^^u » 

^3 =(730-3702)^ +(721 -^^Ozf > =(703 +712^ +(721 +703)^ 
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The feature list including features ' in the iris 
analyzing region is generated based on region segmentation, 
moment invariants are calculated for each feature. The 
moment invariants for effectively discriminating a feature 
5 from another feature exist- Similar images moved, rotated 
and scaled-up/down have similar moment invariants. The 
moment invariants have a difference due to discretization 
error from each other. 

When the size variation of iris is jnodeled as 
10 variation of scale space, if a moment is normalized with a 
mean size, a size-^invariant Zernike moment is generated. 

A radius of the iris image which is transformed to the 
polar coordinates is increased by a predetermined angle, 
the iris image is converted into binary image in order to 
15 obtain a primary contour of the iris having the same radius. 

Histograms are extracted, and it accumulates frequency 
numbers of gray value of pixels in the primary contour of 
the iris in a predetermined angle. In general, to obtain a 
scale space for a discontinuous signal, a continuous 
20 equation should be transformed into a discrete equation by 
using a square formula of integration. 

If F is a smoothen curve of a scale space image, 
wherein the scale space image is scaled by Gaussian kernel, 
a zero-crossing point of a first derivative dF/dK of F in a 
25 scale T is a local minimum value or a local maximum value of 
the smoothen curve in the scale x. A zero-crossing point o^ 
a second derivative d2F/d2x of F is a local minimum value or 
a local maximum value of the first derivative SF/Sx of F in 
the scale x. Ah extreme value of a graSient is_ a point of 
30 inflection in a circular function. The relation between 
the extreme point . and the zero-crossing point is 
illustrated in Fig: 18. 

Referring to Fig. 18, the curve (a) denotes a smoothen 
curve of a scale image in a scale, the function F(x) has 
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three 'extreme points and two minimum points. The 'curve (b) 
denotes zero-crossing points of a first derivate of the 
function F(x) on the extreme points and the minimum points 
of the curve (a) . The zero-crossing points a, c^ e, 
5 indicate the extreme points and the zero-crossing points b, 
d indicate the minimum points. the curve (c) denotes a 
second derivative 52F/52x of the function F and has four 
zero-crossing points f, g, h, . i. The zero-crossing points 
f and h are the minimum values of the first derivate and 
10 starting points of valley regions. The zero-crossing 
points g and i are the extreme values of the first derivate 
and starting points of peak regions. In the range [g, h] , 
a peak region of the circular function is detected. The 
point g is a left gray value and a zero-crossing point of 
15 the second derivate, and the sign of the first derivate 
function on the point g is positive. The point h is a 
right gray value and a zero-crossing point of the second 
derivate, and the sign of the first derivate function on 
the point h is negative. The iris can be represented by 
20 set of the zero-crossing points of the second derivate 
function. Fig. 19 illustrates a peak region and valley 
regions in Fig. 18(a). In Fig. 19, ^^p" denotes a peak 
region, ^^v'' a valley region, a change of sign of the 

'second derivate function from positive to negative, a 
25 • change of sign of the second ^ derivate function from 
negative to positive. A zero contour line can be obtained 
by detecting a peak region ranged from + " to 

According to the above method, an iris curvature 
feature can be illustrated, wherein the iris curvature 
feature represents * shape and movement of inflection points 
of the smoothed signal and is a contour of the zero- 
crossing points of the second derivate. The iris curvature 
feature provides texture of the circular signal in whole 
scales. Based on the iris curvature feature, events 
occurred on the zero-crossing point of a primary contour 
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scale of the shape in the iris analyzing region can be 
detected, the events can be localized by following the- 
zero-crossing points in fine scale step-by-step. A zero 
contour. of the iris curvature feature has a shape of arch, 
5 wherein top portion of the arch is close and bottom portion 
of the arch is open. The zero-crossing points are crossed 
on the peak point of the zero contours as opposite signs, 
which means that the zero-crossing point is not disappeared 
but the scale of the. zero-crossing point is reduced, 

10^ A scale space filtering represents scale of the iris 

• by handling size of filter smoothing the primary contour, 
pixel gray values of the feature in a iris analyzing region 
as a continuous parameter. The filter used for the scale 
space filtering is a filter generated by combining a 

15 Gaussian function and a scale constant. The size of the 
filter used for the scale space filtering is determined 
based on a scale constant, e.g., a standard deviation. The 
size of the filer is expressed as a following equation 40. 

fix. V, O = fix, ^v) ^ six. y,0 

20 Eq. 40 

In the equation 40, V={x(u), y{u), u^ [0,l]), and u is 
a iris image descriptor generated by making the property of 
the iris image as a gray level and binalizihg the iris 
25 image based bn the threshold T. The function f(x, y) is a 
primary contour pixel gray histogram of the iris to be 
analyzed, g(x, y, t) is a .Gaussian function, (x, y, t) is a 
scale space plane. 

In the scale space filtering, the wider region of two- 
30 dimensional image is smoothed as the scale constant t is 
larger. Second deriva'te of F (x, y, i) can be obtained by 
applying V^g (x, y, t) into f(x, y) , which is expressed in 
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10 



15 



a following equation 41. 



Eq, 41 



In the scale space filtering, as the scale constant x. 
is increased, g (x, ' y, t) is increased, and therefore, it 
takes a lot of time to obtain a scale space image- This 
problem can be solved by applying hi, h2, which is expressed 
in a following equation 42- 



X 



The second derivate of F (x, y, t) is expressed in a 
following equation 43. 

: - ■ ' ■ '^l-:- 

In a region in which a result value obtained based on 
V^g (x, y, I) is negative, as the scale space filtering 
20 constant ^,is small, a^ plurality of meaningless peaks are 
generated and the number of the peaks are increased. 
However, if the scale filtering constant is large, e.g., 40, 
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the filter includes the two-dimensional histogram and the 
peak has a shape of combining a plurality of peaks, the 

^ scale space filtering for a larger scale does not effect to. 
find an outstanding peak of the two-dimensional histogram, 
5 In the region in which the values of x and y are larger 
than |3t|, V^g (x, y, x) has a verysmall value . which does 
not affect the calculation result. Therefore, V^g (x, y, 
X) is calculated in a range from -3t to 3t. An image of 
which peak is extracted from second differential of the 

10 scale space image is referred to as a peak image. 

Hereinafter, an automatic optimal scale selection will 
be explained. 

A peak image, which includes outstanding peaks of the 
two-dimensional histogram and represents the shape of the 
15 histogram well, is selected, a scale constant at that time 
is detected at the graph, and then the optimal scale is 
selected. The change of the peak includes four cases as: 
® Generation of a new peak 

® Division of a peak into a plurality of peaks 
20 ® Combination of a plurality of peaks into a new peak 

® Change of shape of peak 



The peak is represented as a node in the graph, and 
relation between peaks of two adjacent peak images^ is 
represented by a directional peak. The node includes a 
scale constant at which the peak starts and a counter, a 
range of scale .in which the peak continuously appears is 
recorded, and a range of scale in which the outstanding 
peaks simultaneously exist is determined. 

A start node is generated, nodes for the peak image 
corresponding to the scale constant 40 are generated, when 
the change of the peak corresponds to the case ®. ® or ®, a 
new node is generated, a start scale of the new node is 
recorded and the counter is operated. If the graph is 
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completed, all of paths from the start node to a 
termination node are searched, a scale range of an 
outstanding peak in each path is founded. In case that a 
new peak is generated, a valley region in the previous peak 
5 image is changed into a peak region due to the change of 
the scale- If there is only one peak newly generated in a 
path and the scale range of the peak is larger than the 
scale range of the valley, since the peak can not t>e 
regarded as an outstanding peak, the scale range of the 
10 outstanding peak is not founded. The range in which the 
scale ranges aire overlapped is determined as a variable 
range, the smallest scale constant within the variable 
range is determined as the optimum scale. (See Fig. 20) 

Hereinafter, a shape descriptor extracting procedure 
15 S313 will be described. 

A shape descriptor extractor 29 generates a Zernike 
moment based on features points extracted from the scale 
space and the scale illumination, and extracts based on the 
Zernike moment a shape descriptor which is rotation- 
20 invariant and strong to an error. At this time, 24 
absolute values of the 8^^ Zernike moment are used as the 
shape descriptor in order to solve the problem that the 
Zernike moment is sensitive to the size of the image and 
the light, by using the scale space and the scale. 
25 The shape descriptor is extracted based on the 

normalized iris curvature ' feature obtained in the pre- 
processing procedure. Since the Zernike moment is 
extracted based on internal region of the iris curvature 
feature, and is rotation-invariant and strong to an error; 
30 the Zernike moment is widely used for a pattern recognition 
system. In this embodiment, as a shape descriptor for 
extracting shape information from the normalized iris 
curvature feature, 24 absolute values of the first to the 
8*"** Zjsrnike moments except of the 0^*^ moment. Also, movement 
35 and scale normalization affect on two Zernike moments Aoo 
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k mOO 



and All. In the normalized "image, there are lAool = (2/ic)^ 

= l/nand I An 1=0. 

Since each of lAool and I An I has the same value in all 
• of the normalized' images, the" moments are excluded from" 
5. feature vector used for representing the features of the. 
image. The 0*"*" moment represents the size of the image and^ 
is used for obtaining a size-invariant feature value. By 
modeling the variation in the size of the image based on 
variation in the spale space, the moment is normalized as a 
10 mean size, to thereby generate the Zernike moment. 

,The Zernike moment f (x, y) of two-dimensional image is 
a complex - moment defined by a following equation 45. The 
Zernike moment is known to have rotation-invariant feature. 
The Zernike moment is defined as a complex polynomial set 
15 each of which element is orthogonal within a unit circle 
(x^+y^^l) . The complex polynomial set is defined as a 
following equation 44. 



20 



zp = lYnm{x,y)\x^ + ^ l) Eq . 4 4 



A basis function of the Zernike moment is. expressed by 
a following equation 45, a rotational axis is a complex 
function defined within a unit circle (x^+y^^l), Rnm(p) is an 
orthogonal radial polynomial equation. Rnin(p) is defined as 
25 the equation 45. 

K»(^.:>') = K.„(A^) = ^»,(P)«""'' Eq. 45 

Where the condition n-|m|: even number and |m|^n 
30 should be satisfied when n is an integer equal to or larger 
than 0, m is an integer. 

In other words, degree n is repeated by m, which is 

expressed as: , p"'^ , p*"* . Wherein /> = V^'+/ , 
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^ = tan"*(— ) . 0 represents an angle between x-axis and the 

X 

vector y. 

Rnm(p) is polar coordinates of Rnm(x^ y) • 
Rnm(p) is a polar coordinate of Rnm(x, y), i.e., x = pcos0, 
5 >' = psiciO . . 

^2 2 

Where, Rn,-m() is equal to Rnxa(p). = p''"'i'/''^^(2p' -1) is 

10 a Jacobi's polynomial equation under a condition that 
s = (n.|m|)/2, P/'J'^Cx). 

A recursive equation of Jacobi's polynomial is used 
for calculating Rnm(p) in order to calculate Zernike 
polynomial without a look-up table. 
15 The Zernike moment for iris curvature feature f(x, y) 

obtained from iris within a predetermined angle by a scale- 
space filter is a Zernike orthogonal basis function, i.e., 
a projection of f(x, y) for Vnni(x, y) . Applying the n-th 
Zernike moment satisfying p", p"*^ / ..••/ p'"' to a discrete 
20 function (not a continuous function) , the Zernike moment is 
a complex number calculated by the Zernike moment expressed 
by an equation 47. 

25 

Wherein * denotes a complex conjugate of \y^{x^yij\. 
A^=—YLf^'''y^^'^nn.ix.y)^jVI„„ix,y)\, x^+y^:^! Eq. 48 
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Wherein VR is a real component of [V^{x^y)X and vi an 

imaginary component of [V^^X^y)] . 

If the Zernike moment for the iris curvature feature 
• 5 f(x, y) is . a Zernike moment (expressed by an equation 
49) of a rotated signal is defined as equations 50 and 51. 

/' (p, 0) = f(P, a + 0) = F(y cos(a) + x sin(a), y sin(cr) - x cos(a)) 

Eq. 49 

10 ^^=— SZ/(Aa + ^)»''"-(P.^). s,t,p^l Eq. 50 

A'^=A^cxp(.-Jma) Eq. 51 

M'*-| = M,J . Eq. 52 

As shown in Eq. 52, an absolute value of the Zernike 
15 moment has the same value regardless rotation of the 
feature. In real computation, if the order of the moments 
is too low, the patterns are difficult to be classified, 
and if the order of the moments is too high, the amount of 
the computation is too large. It is preferable that the 
20 order of the moment is 8 (Refer to Table 4). 

Table 4 
IA„| 

IA20I, IA22I 
lAail.lAaal 

-I-A40I . I A42IVI A44-I ^ — ■ ■ 

IA51L IA53I. IA55I 
lAeol. IA62MA64I.. lAeel 
IA71I. IA73I.IA75I, IA77I 
lAsol. iAaal. lAMl.lAsel. lAssl 



Since the Zernike moment is calculated based on the 
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orthogonal polynomial equation^ the'Zernike moment has a. 
rotation-invariant feature. In particular, the Zernike 
moment has better characteristics in iris representation, 
. duplication and noise. However, the Zernike moment has 

5 * shortcomings to be sensitive to the size and the brightness 
of the image. The shortcoming related to the size of the 
image can be solved based on the scale-space of the image. 
Using Pyramid algorithm, a pattern of the iris is destroyed 
due to the re-sampling of the image. However, the scale- 

,0 space algorithm has better feature point extraction 
characteristic than the Pyramid algorithm, because the 
scale-space algorithm uses the Gaussian function. 
Modifying the Zernike moment, which is invariant to 
movement, rotation and scale of the image, can be extracted 

.5 (refer to an equation 53) . In other words, the image is 
smoothed based on the scale-space algorithm, and the 
smoothed image is normalized, the Zernike moment is robust 
to the size of the image. 



20 = ^ iilos\F,ip,Of ^^^^^0^pdpd0 Eq. . 53 

The modified rotation invariant transform has a 
characteristic that a low , frequency component is emphasized. 
On the other hand, when modeling local luminance variation 

25 -expressed by an equation 54, a JDrightness-invariant Zernike 

moment as expressed by an equation 55 can be generated by 
normalizing the moment by a mean brightness ZOO. 



Eq. 54 



30 
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Z(f'ix,y)) aMx,y)) _ zjf jx.y)) ^^ 55 

Wherein f (x^ y) denotes an iris image^ f^(x^ y) an 
iris image under a new Iviminance, a^ a local luminance 
5 variation rate^ mf a mean Iximinance (a mean, luminance of 
the smoothed image), and Z a Zernike moment operator. 

Though the iris image inputted based on the above 
features is modified by the movement , the scale and the 
rotation of the iris image, the iris pattern, which is 

10 modified in a similar as visual characteristics of the 
human being, can be retrieved. In other words, the shape 
descriptor extractor 29 of Fig, 2 extracts features of the 
iris image from the input image, the reference value 
storing unit 30 of Fig. 2 or the iris pattern registering 

15 unit 14 of Fig. 1 stores the features of the iris image on 
the iris database (DB) 15 at steps S314 and S315. 

If a query image is received at step S316, the shape 
descriptor extractor 2 9 of Fig. 2 or the iris pattern 
feature extractor 13 of Fig. 1 extracts shape descriptors 

20 of the query image (hereinafter, which is referred to as a 
^'query shape descriptor") . The iris pattern recognition 
unit 16 compares the query shape descriptor and the shape 
descriptors stored on the iris DB 15 at step 5317, 
retrieves images corresponding to the shape descriptor 

25 having the minimum distance from the query shape descriptor, 
and outputs the retrieved image to the user. The user can 
see the retrieved iris images rapidly. 

The steps S314, 315 and 317 will . be described in 
detail. 

30 The reference value storing unit 30 of Fig. 2 or the 

iris pattern registering unit 14 of Fig. 1 classifies the 
images as template type based on stability of the Zernike 
moment and similarity according to a Euclidean distance, 
and stores the features of the iris image on the iris 
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database (DB) 15 at step S314, wherein the stability of 
Zernike moments relates to sensitivity which is four- 
directional standard deviation of the Zernike moment. In 
other words, the image patterns of the iris curvature f (x, 

5 y) are projected to the Zernike complex polynomial equation 
Vnm(x, y)on 25 spaces, and classified. The stability is 
obtained by comparing feature points of the current image 
and the previous image, i.e., comparison of the locations 
of the feature points. The similarity is obtained by 

10 comparing distance of areas. Since there are many 
components of the Zernike moment, the area is not a simple 
area, the component is referred to as a template. When 
defining an image analysis region of the image, sample data 
of the image is gathered. Based on the sample data of the 

15 image, the similarity and the stability are obtained. 

The image recognizing/verifying unit 31 of Fig. 2 or 
the iris pattern recognition unit 16 of Fig. 1 recognizes a 
similar iris image by matching features of models which are 
modeled based on the stability and the similarity of the 

20 Zernike moments, and verifies the similar iris image based 
on a least square (LS) algorithm and a least median of 
square (LmedS) algorithm. At this time^ the distance of 
the similarity is calculated based on Minkowsky Mahalanbis 
distance. 

25 The present invention provides a new similarity 

measuring method appropriate for extracting feature 
invariant to the size and luminance of the image, which is 
generated by modifying the Zernike moments.^ 
^ The iris recognition system includes a feature 

30 extracting unit and a feature matching unit. 

In the off-line system, the Zernike moment is 
generated based on the feature point extracted in the scale 
space for the registered iris pattern. In real time 
recognition system, the similar iris pattern is recognized 

35 by statistical matching of the models and the Zernike 
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moment generated based on the feature point, and verifies 
the similar iris pattern by using the LS algorithm and the 
LmedS algorithm. 

Classification of iris images into templates will be 
5 • described in detail. 

In the present invention, the statistical iris 
recognition method recognizes the iris by reflecting the * 
stability of the Zernike moment and the , similarity of 
characteristics to the model statistically. 
10 A basis definition of the modeling is followed. 

An input image is denoted by S, a set of models M={Mi}, 
i=l, 2, Nm, wherein Nm is the number of the models, a set 
of the Zernike moments of the input image S Z={Zi}, i=l, 2, 
.../ Nsf wherein Ns is the number of the Zernike moments of 
15 the input image S- The Zernike moment of the model 
corresponding to the i"~th Zernike moment of the input image 
S is expressed as Zi={Z/ }, j = 1, 2, Nc, wherein Nc is 

the number of the corresponding Zernike moments. 

The probability iris recognition finds a mode Mi which 
20 makes a maximum probability value when the input image S is 
received, which is expressed by an equation 56. 

argmax/*(A^,|5) Eq. 56 

25 A hypothesis as a following equation 57 can be made 

based on candidate model Zernike moments corresponding to 
the Zernike moments of the input image. 

30 Eq. 57 

Where N„ denotes the number of elements of product of 
the model Zernike moments corresponding to the input image. 
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The tol^l hypothesis set can be expressed as: 

H = {H,kjH,-H^J Eq. 58 

5 Since the hypothesis H includes candidates of the 

features extracted from the input image S can be 

replaced by H. If Bayes' theory is applied to the equation 
56, an equation 59 can be obtained as: 

If the probability that each of irises is inputted is 
the same and independent from each other, the equation 59 
can be expressed by an equation 60. 



15 



In Eq. 60, according to theorem of total probability, 
the denominator can be expressed as: P{Hi)^^P{Hi\M^)P{Mt) 

20 In the equation, it is most important to obtain a 

value of the probability ^(7/^ I Af,) . In order to define a 
transcendental probability* /^(/Ta I ^f^) ^ ^ concept on the 

stability is. introduced.. 

The transcendental probability | M,) has a large 

25 value when the stability cr^' and the similarity tzr^ are large 
The stability represents incompleteness of the feature 
points, and the similarity is obtained by the Euclidean 
distance between the features. 

First, the stability will be described in detail. 
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The stability of the Zernike moments is inverse 
proportion to a sensitivity of the Zernil^fe moment to 
variation in the location of the feature points. The 
sensitivity of the Zernike moment represents standard 
deviation of the Zernike moment in four directions from the 
center point- The sensitivity of the Zernike moment is 
expressed by a following equation 61, The stability of the 
Zernike moment is inverse proportion to the sensitivity of 
the Zernike moment. As the sensitivity of the Zernike 
moment is lower^ the stability of location error of the 
feature point is higher. 

SENSITIVITY = ^[|| - Z Jl^ + 1| Z, - Z + H' Z, - Z, ||^ ] 

Eq. 61 

Next, the similarity cr^ will be described in detail. 

As the Euclidean distance from the model feature 
corresponding, to the Zernike moment of the input image is 
shorter, the similarity cr^ is larger. The similarity cr^ is 
expressed by a following equation 62. 

fZ7 oc ^ Eq« 62 

distance 

The recognition result, can be obtained by 
classification of the patterns after performing pre- 
processing, e.g., normalization, which is expressed by a 
following _equat ion. 63_ as:: _ _.. 



A^=—Y.Ef(''MVR^ix.y)^JVI„„^^^^ x^V^l Eq. 63 

If n=0, 1, 8, m=0, 1, 8, the area pattern of the 
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iris curvature f (x, y) is projected to the Zernike complex 
polynomial V^{x^y) on 25 spaces, X=(xi, X2/ Xm) and G=(gi, 
g2^ "9 gm) are classified as a template in the database and 
stored. The distance frequently used for the iris 
5 recognition is classified as a Minkowsky Mahalanbis 
distance. 

£)(^,G)^[X|x,-g, r Eq. 64 

1=1 

0 Where jc,. denotes a magnitude of the i-th Zernike moment 

of the image stored on the DB, and a magnitude of the i- 

th Zernike moment of the query image. 

In case of q=25, the image having the shortest 
Minkowsky's distance within a predetermined permitted limit 

5 is determined as the iris image corresponding to the query 
image. If there is no image having the shortest 
Minkowsky's distance within the predetermined permitted 
limit, it is determined that there is no studied image 
circular shape. For only easy description, it is assumed 

10 that there are two iris images in the dictionary. 
Referring to Fig. 23, input patterns of the iris image, 
wherein the first and the second ZMMs of the rotated iris 
images in a two-dimensional plane, are located on points a 
and b. Euclidean distances da' a, da'b between the points a 

15 and b ar^ obtained^ based . on a following equation 65, 
wherein the Euclidean distance is an absolute distance in 
case of q=l. Euclidean distances are da'a<da'b, da'a<A, 
' which shows that the iris images are rotated. However, if 
the iris images are the same, ZMMs of the iris images are 

to identical with the predetermined permitted limit. 



/>(^.c7)=|;ix,-g, r 
Ml 
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10 



For retrieving the iris image^ shape descriptors of 
the query image and images stored in the iris database 15 
are extracted, and then the iris image similar to the query 
image is retrieved based on the shape descriptors. The 
distance between the query image and the image stored in 
the iris database 15 is obtained based on a following 
equation 66 (Euclidean distance in . case of q=2) , the 
similarity S is obtained by a following equation 67. 

Z)(jr,G)= SCx,-^,)^ Eq. 66 

V <-i 

S = -^ Eq. 67 

l + £> 

The similarity S is normalized and becomes a value 
15 between 0 and 1. Accordingly, the transcendental 

probability P{H,,\M,) can be obtained based on the stability 
and the similarity, which is expressed by a following 
equation 68 . 

20 P(/fJM,) = jfi>((Z,,Z,|M,) Eq. 68 



Where P(iZ^,Zj\M,) is defined as: 



disKZ^Z,) ^ ,:^.Z^eZ(M,) Eq. 69 



e else 

25 Where Ns is the number of interest points of the input 

image, a is a normalization factor obtained by multiplying 
a threshold of the similarity and a threshold of the 
stability, and e is assigned if the corresponding model 
feature does not belong to a certain model. In this 
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embodiment, e is 0.2. To find matching pairs, it is used 
an approximate nearest neighbor (ANN) search algorithm, ' 
which takes log time for linear search space. 

To finid a solution increasing the probability, a 
5 verifying procedure of the retrieved image based on LS and 
. LmedS algorithms will be described. 

The retrieved iris is verified by' matching the input 
image and the model images. Final feature of the iris can 
be obtained through the verification. To find accurate 
10 matching pairs, the image is filtered based on the 
similarity and the stability used for probabilistic iris 
recognition, and outlier is minimized by regional space 
matching. 

Fig. 24 is a diagram showing a method for matching 
15 local regions based on area ratio in accordance with an 
embodiment of the present invention. 

AP P P 

In continuous four points, values .2 J J model 
AP P* P ' 

and ^ ^ ^ for the input image are obtained, if the ratio 
AP\ P\ P\ 

of the two values is larger than the permitted value, the 
20 fourth pair is deleted. At this time, three pairs are 
assumed to be matched. 

Homography is obtained based on the matching pairs. 
The. homography is calculated based on the least square (LS) 
algorithm by using at least three pairs of feature points. 
25. The homography which makes the outlier a minimum value is 
selected as an initial value, and the homography is 
optimized based on the least median of square (LmedS) 
algorithm.' The models are aligned to the input images 
based on the homography. . If the outlier is over 50%, align 
30 of the models is regarded as fail. As the number of 
matched models is larger than the number of the other 
models, the recognition capacity becomes higher. Based on 
this- feature, a discriminative factor is proposed. The 
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disc'rimiriative factor (DF) is defined as: 

• ,DF--^ . Eq,. 70. . 

No 

5 Where Nc is the number of the matching pairs of the 

^models identical to the query iris image, is the number 
the matching pairs of the other models. 

DF is an important factor to select factors of the 
recognition system. The order of the Zernike moments for 
10 the image having the Gaussian noise (of which the standard 
deviation is 5) is 20. When the size of the local image of 
which center point is a feature point is 21X21, the DF has 
the largest value. 

The retrieval performance of the iris recognition 
15 system will be described. 

To evaluate the performance of the iris recognition 
system, a plurality of iris images are necessary. 
Registration and recognition for a certain person are 
necessary, and the number of the necessary iris images is 
20 increased. Also, since it is important the experiment for 
the iris recognition system in various environments in the 
sexual distinction, age and wearing glasses, in order to 
obtain accurate performance result of the recognition 
experiment, a fine plan for the experiments is necessary. 
25 In this embodiment, iris images of 250 persons are 

used, wherein the iris images are captured by a camera. 
500 false acceptance rate (FAR) images for .registering 250 
users (left and right irises) and 300 false rejection rate 
(FRR) images obtained- from 15 users are used in- this 
30 embodiment. However, image acquisition according to time 
and environment should be studied. Table 5 shows data used 
for evaluating the performance of the iris recognition 
system. 
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Table 5 



Number of 
Users 


250" 


Male 


Female 


168 


82 


. Wearing 
Glasses 


Wearing' Contact 
lenses 


Not 
Wearing 


44 


16 


190 


Obtained Data/ 
User 


(FAR: 500) (FRR: 3^0) 


Total Number 
of Data 


250*2=500, 15*20=300 



The pre-processing procedure is very important to 
improve performance of the iris recognition system* 



Table 6 



Procedure 


Fl 


F1+F2 


F1+F2+F3 


Processing Time 


0.1 


0.2 


0.4 



Fl: grids detection 

F2: pupil location detection 

F3: edge component detection 



Table 7 





Male 


Female 


Total 


Not Wearing 


110 


80 


190 


Wearing Glasses 


10 


34 


44 


Wearing Contact 
lenses 


' 8 


' 8 


16 


Total" 


168 


82 


230 



Table 8 





Number 


Rate (%) 


Normal 
Images 


Normal detection 


500 


100 


100 


Inner boundary detection fail 


0 


0 


Outer boundary detection fail 


0 


0 
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Abnormal 
Images 


Shortage of boundary detection 
information _ . 


0 


0 


0 


Error image 


0 


0 


Total ■ . 


$00 


100 



In general, the recognition system is evaluated by two 
, error rates. The two error rates include a false rejection 
rate (FRR) and a false acceptance rate (FAR) . ERR is a 
5 probability in which a user * fails to authenticate 
himself /herself when trying to authenticate by using 
his/her iris images. FAR is a probability in which another 
user success to authenticate himself /herself when trying to 
authenticate by using his/her iris images. In other words, 
10 in order that the biometric recognition system provides a 
highest stability, the biometric recognition system should 
recognize the registered user accurately when the 
registered user tires to be authenticated, and the 
biometric recognition system should deny the unregistered 
15 user when the unregistered user tires to be authenticated. 
These principles of the biometric recognition system should 
be also applied to the iris recognition system. 

According to application field of the iris recognition 
system, the error rates can be selectively adjusted, 
20 However, to increase performance of the iris recognition 
system, both of the two error rates should be decreased. 

A calculating procedure of the error rates will be 
described. 

After calculating distances between the iris images 
25 acquired from the same person based on a similarity 
calculation method, a distribution of frequencies in the 
distances is calculated, which is referred to as 
''authentic". A distribution of frequencies in the distances 
between the iris images acquired from different persons is 
30 calculated, which is referred as to ''imposter''. . Based on 
the authentic and the imposter, a boundary value minimizing 
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the FRR and .the FAR is calculated. The boundary value is 
referred to as ""threshold" . The st.udied. data are used for 
the above procedures. The FRR and the FAR according to "the 
distribution are illustrated in Fig. 25. 

■5 

number of accepted imposter claims ^^qq^^ Eq 71 

r total number of imposter accesses 

number of rejected client claims ^ ^qq^,^ 
total number of client accesses 

The procedure calculating the two error rates for the 

10 iris recognition system will be described. 

If the distance between the studied data and the iris 
image of the same user is smaller than the threshold^ the 
user is authenticated. However, if the distance is larger 
than the threshold, the iris image is determined to be 

15 different from the studied data and the user is denied. 
These procedures are repeated, the number of rejected 
client claims to the total number of client accesses is 
obtained as FRR. 

FAR is calculated by comparing the studied data with 

20 the iris images of the unregistered user. In other words, 
the registered user is compared with another user 
unregistered. If the distance between the studied data and 
the iris image of the user is smaller than the threshold/ 
the user is determined as the same person. However, if the 

25 distance is larger than the threshold, the . user is 
determined as a different person.' These procedures are 
repeated, the number of acciepted imposter claims, to the 
total number of imposter accesses is obtained as FAR. ^ 

In the present invention, for the verification 

30 performance evaluation, FAR and FRR are performed on data 
selected at the pre-processing. 

The authentic distribution and the imposter 
distribution will be described. 
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After calculating distances between the iris images 
'Acquired from the same person based on a similarity 
calculation method," a distribution of frequencies in 
distances is calculated^ which * is referred to as 
5 "'authentic". The authentic distribution is illustrated in 
Fig, 26. In this drawing, an x-axis. denotes a distance and 
a y-axis a frequency. 

Fig. 27 is a graph showing a distribution in distances 
between iris images of different persons where an x-axis 
10 denotes a distance and a y-axis a frequency. 

It will be described selection of thresholds for the 
authentic distribution and the imposter distribution. 

In general, FRR and FAR are varied according to the 
threshold and can be adjusted according to the application 
15 field. The threshold should be carefully adjusted. 

Fig. 28 is a graph showing an authentic distribution 
and an imposter distribution. 

The threshold is selected based on the authentic 
distribution and the imposter distribution. The iris 
20 recognition system performs authentication based on the 
threshold of an equal error rate (EER) . The threshold of 
the EER is calculated by a following equation 72 expressed 
as : 

25 Threshold = -^^—^^i '—^ Eq. 72 

<T^: standard deviation of authentic distribution 
CTi- sj:andard deviation of Imposter distribution . • 
//^:mean of authentic . distribution 
30 //, :mean of Imposter distribution 

the iris data are classified into studied data and 
text data, and the experiment result is represented in 
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Table 9. • . • . 'f^^ 

^ It takes about 5 to 6 seconds for registrati^if' of the 

image and about 1 to 2 seconds for authent^catipn of the 
query image. 



Table 9 



FRR 


5% 


FAR 


15% 



The present invention can be implemented and stored in 
a computer readable recording medium, e.g., CD-ROM, a 
random access memory (RAM) , a read only memory (ROM) , a 
floppy disk, a hard disk, and a magneto-optical disk. 

While the present invention has been described with 
respect to certain preferred embodiments, it will be 
apparent to those skilled in the art that various changes 
and modifications may be made without departing from the 
scope of the invention as defined in the following claims. 
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